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depth, r, and a slope for the scalar perturbation spectrum, n s ) fits not only 
the three year WMAP temperature and polarization data, but also small scale 
CMB data, light element abundances, large-scale structure observations, and the 
supernova luminosity/distance relationship. Using WMAP data only, the best 
fit values for cosmological parameters for the power-law flat ACDM model are 
(tt m h 2 ,tt b h 2 ,h,n s ,T,a 8 ) = 0.1277^079, 0.02229 ± 0.00073, 0.7321°;°^, 0.958 ± 
0.016,0.089 ± 0.030,0.7611^). The three year data dramatically shrinks the 
allowed volume in this six dimensional parameter space. 

Assuming that the primordial fluctuations are adiabatic with a power law 
spectrum, the WMAP data alone require dark matter, and favor a spectral in- 
dex that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant 
spectrum (n s = l,r = 0). Adding additional data sets improves the con- 
straints on these components and the spectral slope. For power-law models, 
WMAP data alone puts an improved upper limit on the tensor to scalar ra- 
tio, r .oo2 < 0.65 (95% CL) and the combination of WMAP and the lensing- 
normalized SDSS galaxy survey implies r 00 02 < 0.30 (95% CL). 

Models that suppress large-scale power through a running spectral index or a 
large-scale cut-off in the power spectrum are a better fit to the WMAP and small 
scale CMB data than the power-law ACDM model; however, the improvement in 
the fit to the WMAP data is only A% 2 = 3 for 1 extra degree of freedom. Models 
with a running-spectral index are consistent with a higher amplitude of gravity 
waves. 

In a flat universe, the combination of WMAP and the Supernova Legacy 
Survey (SNLS) data yields a significant constraint on the equation of state of the 
dark energy, w = —0-9Q7tomt- If we assume w — — 1, then the deviations from 
the critical density, Q K , are small: the combination of WMAP and the SNLS 
data imply Q k = —0.011 ± 0.012. The combination of WMAP three year data 
plus the HST key project constraint on Hq implies Q k = —0.014 ± 0.017 and 
Q\ = 0.716 ± 0.055. Even if we do not include the prior that the universe is flat, 
by combining WMAP, large-scale structure and supernova data, we can still put 
a strong constraint on the dark energy equation of state, w = —1.08 ± 0.12. 

For a flat universe, the combination of WMAP and other astronomical data 
yield a constraint on the sum of the neutrino masses, < 0.66 eV (95% CL). 

Consistent with the predictions of simple inflationary theories, we detect no signif- 
icant deviations from Gaussianity in the CMB maps using Minkowski functionals, 
the bispectrum, trispectrum, and a new statistic designed to detect large-scale 
anisotropics in the fluctuations. 
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Subject headings: cosmic microwave background, cosmology: observations 

1. Introduction 

The power-law ACDM model fits not only the Wilkinson Microwave Anisotropy Probe 
(WMAP) first year data, but also a wide range of astronomical data (Bennett et al. 2003; 
Spergel et al. 2003). In this model, the universe is spatially flat, homogeneous and isotropic 
on large scales. It is composed of ordinary matter, radiation, and dark matter and has a 
cosmological constant. The primordial fluctuations in this model are adiabatic, nearly scale- 
invariant Gaussian random fluctuations (Komatsu et al. 2003). Six cosmological parameters 
(the density of matter, the density of atoms, the expansion rate of the universe, the amplitude 
of the primordial fluctuations, their scale dependence and the optical depth of the universe) 
are enough to predict not only the statistical properties of the microwave sky, measured by 
WMAP at several hundred thousand points on the sky, but also the large-scale distribution 
of matter and galaxies, mapped by the Sloan Digital Sky Survey (SDSS) and the 2dF Galaxy 
Redshift Survey (2dFGRS). 

With three years of integration, improved beam models, better understanding of sys- 
tematic errors (Jarosik et al. 2006), temperature data (Hinshaw et al. 2006), and polarization 
data (Page et al. 2006), the WMAP data has significantly improved. There have also been 
significant improvements in other astronomical data sets: analysis of galaxy clustering in the 
SDSS (Tegmark et al. 2004a; Eisenstein et al. 2005) and the completion of the 2dFGRS (Cole 
et al. 2005); improvements in small-scale CMB measurements (Kuo et al. 2004; Readhead 
et al. 2004a,b; Grainge et al. 2003; Leitch et al. 2005; Piacentini et al. 2005; Montroy et al. 
2005; O'Dwyer et al. 2005), much larger samples of high redshift supernova (Riess et al. 2004; 
Astier et al. 2005; Nobili et al. 2005; Clocchiatti et al. 2005; Krisciunas et al. 2005); and 
significant improvements in the lensing data (Refregier 2003; Heymans et al. 2005; Semboloni 
et al. 2005; Hoekstra et al. 2005). 

In §2, we describe the basic analysis methodology used, with an emphasis on changes 
since the first year. In §3, we fit the ACDM model to the WMAP temperature and polariza- 
tion data. With its basic parameters fixed at z ~ 1100, this model predicts the properties 
of the low redshift universe: the galaxy power spectrum, the gravitational lensing power 
spectrum, the Hubble constant, and the luminosity-distance relationship. In §4, we compare 
the predictions of this model to a host of astronomical observations. We then discuss the 
results of combined analysis of WMAP data, other astronomical data, and other CMB data 
sets. In §5, we use the WMAP data to constrain the shape of the power spectrum. In §6, 
we consider the implications of the WMAP data for our understanding of inflation. In §7, 
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we use these data sets to constrain the composition of the universe: the equation of state of 
the dark energy, the neutrino masses and the effective number of neutrino species. In §8, we 
search for non-Gaussian features in the microwave background data. The conclusions of our 
analysis are described in §9. 



The basic approach of this paper is similar to that of the first-year WMAP analysis: 
our goal is to find the simplest model that fits the CMB and large-scale structure data. 
Unless explicitly noted in §2.1, we use the methodology described in Verde et al. (2003) and 
applied in Spergel et al. (2003). We use Bayesian statistical techniques to explore the shape 
of the likelihood function, we use Monte Carlo Markov chain methods to explore the likeli- 
hood surface and we quote both our maximum likelihood parameters and the marginalized 
expectation value for each parameter in a given model: 



where al is the value of the i— th parameter in the chain and j indexes the chain element. 
The number of elements (M) in the typical merged Markov chain is at least 50,000 and is 
always long enough to satisfy the Gelman & Rubin (1992) convergence test with R < 1.1. 
In addition, we use the spectral convergence test described in Dunkley et al. (2005) to 
confirm convergence. Most merged chains have over 100,000 elements. We use a uniform 
prior on cosmological parameters, p(a) unless otherwise specified. We refer to < «j > as the 
best fit value for the parameter and the peak of the likelihood function as the best fit model. 

The Markov chain outputs and the marginalized values of the cosmological parameters 
listed in Table 1 are available at http:/ /lambda.gsfc. nasa.gov for all of the models discussed 
in the paper. 



We now use not only the measurements of the temperature power spectrum (TT) and the 
temperature polarization power spectrum (TE), but also measurements of the polarization 
power spectra (EE) and (BB). 

At the lowest multipoles, a number of the approximations used in the first year analysis 
were suboptimal. Efstathiou (2004) notes that a maximum likelihood analysis is significantly 



2. Methodology 




(1) 



2.1. Changes in analysis techniques 



- 5 - 



better than a quadratic estimator analysis at £ = 2. Slosar et al. (2004) note that the 
shape of the likelihood function at £ = 2 is not well approximated by the fitting function 
used in the first year analysis (Verde et al. 2003). More accurate treatments of the low £ 
likelihoods decrease the significance of the evidence for a running spectral index (Efstathiou 
2004; Slosar et al. 2004; O'Dwyer et al. 2004). Hinshaw et al. (2006) and Page et al. (2006) 
describe our approach to addressing this concern: for low multipoles, we explicitly compute 
the likelihood function for the WMAP temperature and polarization maps. For the analysis 
of the polarization maps, we use the resolution N side = 8 iV _1 matrices. This pixel-based 
method is used for Cj T for 2 < £ < 30 and polarization for 2 < £ < 23. For most of the 
analyses in the paper, we use a N side = 8 version of the temperature map for the analysis 
of the low / likelihood that uses a pixel-based version for £ < 12. For the WMAP ACDM 
only case, we use the more time-consuming N side = 16 version of the code. Hinshaw et al. 
(2006) compares various approaches towards computing the low / likelihood. In Appendix 
A, we discuss various choices made in the maximum likelihood code. For the ACDM model, 
we have computed the best fit parameters using a range of assumptions for the amplitude 
of point source contamination and different treatments of the low / likelihood. 

There are several improvements in our analysis of high £ temperature data (Hinshaw 
et al. 2006): better beam models, improved foreground models, and the use of maps with 
smaller pixels (N side = 1024). The improved foreground model is significant at £ < 200. The 
N S ide = 1024 maps significantly reduce the effects of sub-pixel CMB fluctuations and other 
pixelization effects. We found that N side = 512 maps had higher x 2 than N side = 1024 maps, 
particularly for £ = 600 — 700, where there is significant signal-to-noise and pixelization 
effects are significant. 

We now marginalize over the amplitude of Sunyaev-Zel'dovich (SZ) fluctuations. The 
expected level of SZ fluctuations (Refregier et al. 2000; Komatsu & Seljak 2001; Bond et al. 
2005) is £{£ + 1)Q/(2tt) = 19 ± 3^K 2 at £ = 450 - 800 for VL m = 0.26, tt b = 0.044, 
h = 0.72, n s = 0.97 and erg = 0.80. The amplitude of SZ fluctuations is very sensitive 
to og (Komatsu & Kitayama 1999; Komatsu & Seljak 2001). For example at 60 GHz, 
£{£ + l)C e /(2n) = 65 ± 15(/iK) 2 at £ = 450 - 800 for a 8 = 0.91, which is comparable to the 
WMAP statistical errors at the same multipole range. Since the WMAP spectral coverage is 
not sufficient to be able to distinguish CMB fluctuations from SZ fluctuations (see discussion 
in Hinshaw et al. (2006)), we marginalize over its amplitude using the Komatsu & Seljak 
(2002) analytical model for the shape of the SZ fluctuations. We impose the prior that the 
SZ signal is between and 2 times the Komatsu & Seljak (2002) value. Consistent with the 
analysis of Huffenberger et al. (2004), we find that the SZ contribution is not a significant 
contaminant to the CMB signal on the scales probed by the WMAP experiment. We report 
the amplitude of the SZ signal normalized to the Komatsu & Seljak (2002) predictions for 
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the cosmological parameters listed above with cr 8 = 0.80. A sz = 1 implies that the SZ 
contribution is 8.4, 18.7 and 25.2 (/iK) 2 at £ = 220,600 and 1000 respectively. We discuss 
the effects of this marginalization in Appendix A. We have checked that gravitational lensing 
of the microwave background, the next most significant secondary effect after the thermal 
SZ effect (Lewis & Challinor 2006) does not have a significant effect on parameters. 

We now use the CAMB code (Lewis et al. 2000) for our analysis of the WMAP power 
spectrum. The CAMB code is derived from CMBFAST (Zaldarriaga & Seljak 2000), but has 
the advantage of running a factor of 2 faster on the Silicon Graphics, Inc. (SGI) machines 
used for the analysis in this paper. For the multipole range probed by WMAP, the numerical 
uncertainties and physical uncertainties in theoretical calculations of multipoles are about 
1 part in 10 3 (Seljak et al. 2003), significantly smaller than the experimental uncertainties. 
When we compare the results to large-scale structure and lensing calculations, the analytical 
treatments of the growth of structure in the non-linear regime is accurate to better than 10% 
on the smallest scales considered in this paper (Smith et al. 2003). 

2.2. Parameter choices 

We consider constraints on the hot big bang cosmological scenario with Gaussian, adi- 
abatic primordial fluctuations as would arise from single field, slow-roll inflation. We do not 
consider the influence of isocurvature modes nor the subsequent production of fluctuations 
from topological defects or unstable particle decay. 

We parameterize our cosmological model in terms of 15 parameters: 

p = {coi,, uj c , r, Q\, w, fife, f„, N v , A^, n s , r, dn s /dlak, Asz, bsDSS, z s) (2) 

where these parameters are defined in Table 1. For the basic power-law ACDM model, 
we use u>b, uj c , exp(— 2r), O s , n s , and C|^ 20 , as the cosmological parameters in the chain, 
Asz as a nuisance parameter, and unless otherwise noted, we assume a flat prior on these 
parameters. Note that r is the optical depth since reionization. Prior to reionization, x e 
is set to the standard value for the residual ionization computed in RECFAST (Seager, 
Sasselov & Scott 2000). For other models, we use these same basic seven parameters plus 
the additional parameters noted in the text. Other standard cosmological parameters (also 
defined in Table 1), such as cr 8 and h, are functions of these six parameters. Appendix A 
discusses the dependence of results on the choice of priors. 

With only one year of WMAP data, there were significant degeneracies even in the 
ACDM model: there was a long degenerate valley in n s — r space and there was also a 
significant degeneracy between n s and Vl b h 2 (see Figure 5 in Spergel et al. (2003)). With the 
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measurements of the rise to the third peak (Hinshaw et al. 2006) and the EE power spectrum 
(Page et al. 2006), these degeneracies are now broken (see §3). However, more general 
models, most notably those with non-flat cosmologies and with richer dark energy or matter 
content, have strong parameter degeneracies. For models with adiabatic fluctuations, the 
WMAP data constrains the ratio of the matter density /radiation density, effectively, Q m h 2 , 
the baryon density, Qbh 2 , the slope of the primordial power spectrum and the distance to the 
surface of last scatter. In a flat vacuum energy dominated universe, this distance is a function 
only of f2 and h, so the matter density and Hubble constant are well constrained. On the 
other hand, in non-flat models, there is a degeneracy between Q m , h and the curvature (see 
§7.3). Similarly, in models with more complicated dark energy properties (w ^ —1), there 
is a degeneracy between fl m , h and w. In models where the number of neutrino species is 
not fixed, the energy density in radiation is no longer known so that the WMAP data only 
constrains a combination of Q m h 2 and the number of neutrino species. These degeneracies 
slow convergence as the Markov chains need to explore degenerate valleys in the likelihood 
surface. 



3. ACDM Model: Does it still fit the data? 
3.1. WMAP only 

The ACDM model is still an excellent fit to the WMAP data. With longer integration 
times and smaller pixels, the errors in the high / temperature multipoles have shrunk by 
more than a factor of three. As the data have improved, the likelihood function remains 
peaked around the maximum likelihood peak of the first year WMAP value. With longer 
integration, the most discrepant high £ points from the year-one data are now much closer 
to the best fit model (see Figure 2). For the first year WMAP TT and TE data (Spergel 
et al. 2003), the reduced xtff was 1-09 for 893 degrees of freedom (D.O.F.) for the TT data 
and was 1.066 for the combined TT and TE data (893+449=1342 D.O.F.). For the three 
year data, which has much smaller errors for £ > 350, the reduced xtff f° r 982 D.O.F. 
{£ = 13 - 1000- 7 parameters) is now 1.068 for the TT data and 1.041 for the combined TT 
and TE data ( 1410 D.O.F., including TE £ = 24 - 450), where the TE data contribution 
is evaluated from £ = 24 — 500. For the T, Q, and U maps using the pixel based likelihood 
we obtain a reduced xtff = 0-981 for 1838 pixels (corresponding to Cf T for £ = 2 — 12 and 
Cj E for £ = 2 — 23). The combined reduced xtff = 1-037 for 3162 degrees of freedom for 
the combined fit to the TT and TE power spectrum at high £ and the T, Q and U maps at 
low £. 

While many of the maximum likelihood parameter values (Table 2, columns 3 and 7 
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Parameter 


Description 


Definition 


Ho 


Hubble expansion factor 


H = 100ft Mpc^km s" 1 




Baryon density 


uj b = fl b h 2 = p(,/1.88 x 10~ 26 kg m- 3 




Cold dark matter density 


uj c = fl c h 2 = p c /18.8 yoctograms/ m~ 3 


u 


Massive neutrino fraction 


fv — 1 




Total neutrino mass (eV) 


£ to„ = 94fi„ft 2 


N v 


Effective number of relativistic neutrino species 






Spatial curvature 




^DE 


Dark energy density 


For w = — 1, f^A = Qde 




Matter energy density 




W 


Dark energy equation of state 


w = pde/pde 


A 2 


Amplitude of curvature perturbations 1Z 


A^(fc = 0.002/Mpc) w 29.5 x 10~ 10 A 


A 


Amplitude of density fluctuations (k = 0.002/Mpc) 


See Spergel et al. (2003) 


n s 


Scalar spectral index at 0.002/Mpc 




a 


Running in scalar spectral index 


a = dn s /dlnk (assume constant) 


r 


Ratio of the amplitude of tensor fluctuations 
to scalar potential fluctuations at k=0.002/Mpc 




TH 


Tensor spectral index 


Assume n t = — r/8 


T 


Reionization optical depth 






Linear theory amplitude of matter 
fluctuations on 8 ft -1 Mpc 




e s 


Acoustic peak scale (degrees) 


see Kosowsky et al. (2002) 


A sz 


SZ marginalization factor 


see appendix A 


bsdss 


Galaxy bias factor for SDSS sample 


b = [P sdss (k, z = 0)/P(fc)] 1 / 2 (constant) 


r< tt 

L/ 220 


Amplitude of the TT temperature power spectrum at I 


= 220 


Zs 


Weak lensing source redshift 





Table 1: Cosmological parameters used in the analysis, http://lambda.gsfc.nasa.gov lists the marginalized 
values for these parameters for all of the models discussed in this paper. 
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Fig. 1. — The improvement in parameter constraints for the power-law ACDM 
model (Model M5 in Table 3). The contours show the 68% and 95% joint 2-d 
marginalized contours for the (Q m h 2 ,a$) plane (left) and the (n s ,r) plane (right). 
The black contours are for the first year WMAP data (with no prior on r). The 
red contours are for the first year WMAP data combined with CBI and ACBAR 
(WMAPext in Spergel et al. (2003)). The blue contours are for the three year 
WMAP data only with the SZ contribution set to to maintain consistency 
with the first year analysis. The WMAP measurements of EE power spectrum 
provide a strong constraint on the value of r. The models with no reionization 
(r = 0) or a scale-invariant spectrum (n s = 1) are both disfavored at ^-xlff > 6 for 
5 parameters (see Table 3). Improvements in the measurement of the amplitude 
of the third peak yield better constraints on Vt m h 2 . 
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200 400 600 800 1000 

Multipole moment I 

Fig. 2. — Comparison of the predictions of the different best fit models to the 
data. The black line is the angular power spectrum predicted for the best fit 
three-year WMAP only ACDM model. The red line is the best fit to the 1-year 
WMAP data. The orange line is the best fit to the combination of the 1-year 
WMAP data, CBI and ACBAR (WMAPext in Spergel et al. (2003)). The solid 
data points are for the 3 year data and the light gray data points are for the first 
year data. 
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and Figure 1) have not changed significantly, there has been a noticeable reduction in the 
marginalized value for the optical depth, r, and a shift in the best fit value of fl m h 2 . (Each 
shift is slightly larger than la). The addition of the EE data now eliminates a large region of 
parameter space with large r and n s that was consistent with the first year data. With only 
the first year data set, the likelihood surface was very flat. It covered a ridge in r — n s over 
a region that extended from r ~ 0.07 to nearly r = 0.3. If the optical depth of the universe 
were as large as r = 0.3 (a value consistent with the first year data), then the measured 
EE signal would have been 10 times larger than the value reported in Page et al. (2006). 
On the other hand, an optical depth of r = 0.05 would produce one quarter of the detected 
EE signal. As discussed in Page et al. (2006), the reionization signal is now based primarily 
on the EE signal rather than the TE signal. See figure 26 in Page et al. (2006) for the 
likelihood plot for r: note that the form of this likelihood function is relatively insensitive 
to the cosmological model (over the range considered in this paper). 

There has also been a significant reduction in the uncertainties in the matter density, 
VL m h 2 . With the first year of WMAP data, the third peak was poorly constrained (see the 
light gray data points in Figure 2). With three years of integration, the WMAP data better 
constrain the height of the third peak: WMAP is now cosmic variance limited up to I = 400 
and the signal-to- noise ratio exceeds unity up to I = 850. The new best fit WMAP-only 
model is close to the WMAP (first year)+CBI+ACBAR model in the third peak region. As 
a result, the preferred value of VL m h 2 now shifts closer to the "WMAPext" value reported 
in Spergel et al. (2003). Figure 1 shows the Vt m h 2 — a 8 likelihood surfaces for the first year 
WMAP data, the first year WMAPext data and the three year WMAP data. The accurately 
determined peak position constrains Q^ 75 h (Page et al. 2003a), fixes the cosmological age, 
and determines the direction of the degeneracy surface. With 1 year data, the best fit value 
is n° m 275 h = 0.498. With three years of data, the best fit shifts to 0.492^;^. The lower 
third peak implies a smaller value of Q m h 2 and because of the peak constraint, a lower value 
of Q m . 

The best fit value of as (marginalized over the other parameters) is now noticeably 
smaller for the three year data alone, 0.761^QQ4g than for the first-year WMAP data alone, 
0.92 ± 0.10. This lower value is due to a smaller third peak height, which leads to a lower 
value of f2 m , and less structure growth and a lower best fit value for r. The height of the 
third peak was very uncertain with the first-year data alone. In Spergel et al. (2003), we 
used external CMB data sets to constrain the third peak and with these data, the Maximum 
Likelihood value was 0.84. With three years of data, the third peak is better determined 
and its height is close to the value estimated from the ground based data. With the EE 
measurements, we have eliminated most of the high r region of parameter space. Since 
higher values of r imply a higher amplitude of primordial fluctuations, the best fit value of 
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a is proportional to exp(r). For a model with all other parameters fixed, r = 0.10 implies a 
7% lower value of cr 8 . As discussed in §4.1, this lower value of a 8 is more consistent with the 
X-ray measurements but lower than the best fit value from recent lensing surveys. Lower a 8 
implies later growth of structure. 

In the first year data, we assumed that the SZ contribution to the WMAP data was 
negligible. Appendix A discusses the change in priors and the change in the SZ treatment 
and their effects on parameters: marginalizing over SZ most significantly shifts n s and o" 8 
by 1% and 3% respectively. In Table 2 in the column labeled No SZ and Figure 1, we 
assume A$z = and use the first year likelihood code to make a consistent comparison 
between the first-year and three- year results. As in the first year analysis, we use a flat prior 
on the logarithm of the amplitude and a flat prior on 6 S and a flat prior on r. The first 
column of Table 5 list the parameters fit to the WMAP three- year data with A sz allowed to 
vary between and 2. In the tables, the "mean" value is calculated according to equation 
(1) and the "Maximum Likelihood (ML)" value is the value at the peak of the likelihood 
function. In the last two columns, we provide our current best estimate of parameters 
including SZ marginalization and using the full N S id e = 16 likelihood code to compute the 
TT likelihoods. In subsequent tables and figures, we will allow the SZ contribution to vary 
and quote the appropriate marginalized values. Allowing for an SZ contribution lowers the 
best fit primordial contribution at high £, thus, the best fit models with an SZ contribution 
have lower n s and os values. However, in other tables, we use the faster N si d e = 8 likelihood 
code unless specifically noted. In all of the Tables, we quote the 68% confidence intervals on 
parameters and the 95% confidence limits on bounded parameters. 



3.2. Reionization History 

Since the Kogut et al. (2003) detection of r, the physics of reionization has been a 
subject of extensive theoretical study (Cen 2003; Ciardi et al. 2003; Haiman & Holder 2003; 
Madau et al. 2004; Oh & Haiman 2003; Venkatesan et al. 2003; Ricotti & Ostriker 2004; 
Sokasian et al. 2004; Somerville & Livio 2003; Wyithe & Loeb 2003; Iliev et al. 2005). Page 
et al. (2006) provides a detailed discussion of the new polarization data: while the best fit 
value for r has not changed significantly, the new EE data, combined with an improved 
treatment of the TE data, implies smaller marginalized maximum likelihood value. The 
three year data favors r ~ 0.1, consistent with the predictions of a number of simulations 
of ACDM models. For example, Ciardi et al. (2003) ACDM simulations predict r = 0.104 
for parameters consistent with the WMAP primordial power spectrum. Tumlinson et al. 
(2004) use the nucleosynthetic data to derive and construct an initial mass function (IMF) 
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Table 2: Power Law ACDM Model Parameters and 68% Confidence Intervals. The Three 
Year fits in the columns labeled "No SZ" use the likelihood formalism of the first year paper 
and assume no SZ contribution, Asz = 0, to allow direct comparison with the First Year 
results. Fits that include SZ marginalization are given in the last two columns of the upper 
and lower tables and represent our best estimate of these parameters. The last column 
includes all data sets. 



Parameter 


First Year WMAPext 
Mean Mean 


Three Year Three Year 
Mean (No SZ) Mean 


Three Year+ALL 
Mean 


ioon b h 2 

Ho 

T 

n s 


9 Qo+0.13 9 09+O.I2 
z - oo -0.12 z - oz -0.11 

n 144+0.016 q 1 o 4 +o.oo6 

72+J 73+j 
n 1 7 +0 - 08 n 1 5 +a07 

u - 1 '-0.07 u - io -0.07 

n qq+ - 04 n 98 +a03 

u - yy -0.04 u - yo -0.03 

29+ - 07 n 2 c i +0 - 03 

u - zy -0.07 u - ZJ -0.03 
u - yz -0.1 u -°^-0.06 


2.23 ±0.08 2.229 ±0.073 
0.126 ±0.009 0.1277±g;ggfg 

73.5 ± 3.2 73.2+Ji 

0.088lg;^ 0.089 ± 0.030 
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0.97 0.96 
0.32 0.27 
0.88 0.82 


2.23 2.22 
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0.0904 0.091 

0.95 0.954 
0.232 0.236 
0.737 0.756 


2.19 
0.131 

73.2 
0.0867 
0.949 
0.259 
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Fig. 3. — WMAP constraints on the reionization history. (Left) The 68% and 95% 
joint 2-d marginalized confidence level contours for x° e — z re i on for a power law 
A Cold Dark Matter (ACDM) model with the reionization history described by 
equation 3 and fit to the WMAP three year data. In equation 3 we assume that 
the universe was partially reionized at 2 re ion to an ionization fraction of x®, and 
then became fully ionized at z = 7. (Right) The 68% and 95% joint 2-d marginal- 
ized confidence level contours for x\ — n s . This figure shows that x® and n s are 
nearly independent for a given value of r, indicating that WMAP determinations 
of cosmological parameters are not affected by details of the reionization history. 
Note that we assume a uniform prior on z re i on in this calculation, which favors 
models with lower x® values in the right panel. 
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for reionization and find r ~ 0.1. Chiu, Fan & Ostriker (2003) found that their joint analysis 
of the WMAP and SDSS quasar data favored a model with r = 0.11, a 8 — 0.83 and n = 0.96, 
very close to our new best fit values. Wyithe & Cen (2006) predict that if the product of star 
formation efficiency and escape fraction for Pop-Ill stars is comparable to that for Pop-II 
stars, r = 0.09 — 0.12 with reionization histories characterized by an extended ionization 
plateau from z = 7 — 12. They argue that this result holds regardless of the redshift where 
the intergalactic medium (IGM) becomes enriched with metals. 

Measurements of the EE and TE power spectrum are a powerful probe of early star 
formation and an important complement to other astronomical measurements. Observations 
of galaxies (Malhotra & Rhoads 2004), quasars (Fan et al. 2005) and gamma ray bursts 
(Totani et al. 2005) imply that the universe was mostly ionized by z = 6. The detection 
of large-scale TE and EE signal (Page et al. 2006) implies that the universe was mostly 
reionized at even higher redshift. CMB observations have the potential to constrain some 
of the details of reionization, as the shape of the CMB EE power spectrum is sensitive 
to reionization history (Kaplinghat et al. 2003; Hu & Holder 2003). Here, we explore the 
ability of the current EE data to constrain reionization by postulating a two stage process 
as a toy model. During the first stage, the universe is partially reionized at redshift z re i on 
and complete reionization occurs at z = 7: 

z > Z re i on 

z<7 (3) 

We have modified CAMB to include this reionization history. 

Figure 3 shows the likelihood surface for x° e and z reion . The plot shows that the data 
does not yet constrain x\ and that the characteristic redshift of reionization is sensitive to 
our assumptions about reionization. If we assume that the universe is fully reionized, x° e — 1, 
then the maximum likelihood peak is z re i on = 11.3. The maximum likelihood peak value of 
the cosmic age at the reionization epoch is t re i on = 365Myr. 

Reionization alters the TT power spectrum by suppressing fluctuations on scales smaller 
than the horizon size at the epoch of reionization. Without strong constraints from polar- 
ization data on r, there is a strong degeneracy between spectral index and r in likelihood 
fits (Spergel et al. 2003). The polarization measurements now strongly constrain r; however, 
there is still significant uncertainty in x e and the details of the reionization history. For- 
tunately, the temperature power spectrum mostly depends on the amplitude of the optical 
depth signal, r, so that the other fit parameters (e.g., n s ) are insensitive to the details of 
the reionization history (see Figure 3). Because of this weak correlation, we will assume a 
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simple reionization history (x° e = 1) in all of the other analysis in this paper. Allowing for a 
more complex history is not likely to alter any of the conclusions of the other sections. 

3.3. How Many Parameters Do We Need to Fit the WMAP Data? 

In this subsection, we compare the power-law ACDM to other cosmological models. We 
consider both simpler models with fewer parameters and models with additional physics, 
characterized by additional parameters. We quantify the relative goodness of fit of the 
models, 

Ax 2 eff = ~ A (2 In £) = 2 In £(ACDM) - 2 In £(model) (4) 

A positive value for Axlff implies the model is disfavored. A negative value means that the 
model is a better fit. We also characterize each model by the number of free parameters, 
N par . There are 3162 degrees of freedom in the combination of T, Q, and U maps and high 
t TT and TE power spectra used in the fits and 1448 independent CVs, so that the effective 
number of data degrees of freedom is between 1448 and 3162. 

Table 3 shows that the power-law ACDM is a significantly better fit than the simpler 
models. For consistency, all of the models are computed with the N sia - e = 8 likelihood code 

Table 3: Goodness of Fit, Axlff = — 21n£, for WMAP data only relative to a Power-Law 
ACDM model. Axlff > is a worse fit to the data. 





Model 


-A(21n£) 


N 

1 y par 


Ml 


Scale Invariant Fluctuations (n s = 1) 


6 


5 


M2 


No Reionization (r = 0) 


7.4 


5 


M3 


No Dark Matter (fi c = 0, fi A ^ 0) 


248 


6 


M4 


No Cosmological Constant (Q c ^ 0, Q\ = 0) 





6 


M5 


Power Law ACDM 





6 


M6 


Quintessence (w ^ —1) 





7 


M7 


Massive Neutrino (m„ > 0) 


-1 


7 


M8 


Tensor Modes (r > 0) 





7 


M9 


Running Spectral Index (dn s /d\nk ^ 0) 


-4 


7 


M10 


Non-flat Universe (Q k ^ 0) 


-2 


7 


Mil 


Running Spectral Index & Tensor Modes 


-4 


8 


M12 


Sharp cutoff 


-1 


7 


M13 


Binned A^(fc) 


-22 


20 
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and the higher value of the point source amplitude. If we reduce the number of parameters 
in the model, the cosmological fits significantly worsen: 

• Cold dark matter serves as a significant forcing term that changes the acoustic peak 
structure. Alternative gravity models (e.g., MOND), and all baryons-only models, lack 
this forcing term so they predict a much lower third peak than is observed by WMAP 
and small scale CMB experiments (McGaugh 2004; Skordis et al. 2006). Models with- 
out dark matter (even if we allow for a cosmological constant) are very poor fits to the 
data. 

• Positively curved models without a cosmological constant are consistent with the 
WMAP data alone: a model with the same six parameters and the prior that there is 
no dark energy, fi A = 0, fits as well as the standard model with the flat universe prior, 
Q m + Q\ — 1. However, if we imposed a prior that H > 40 km s _1 Mpc -1 , then the 
WMAP data would not be consistent with £7a = 0. Moreover, the parameters fit to the 
no-cosmological-constant model, (H = 30 km s _1 Mpc -1 and tt m = 1.3) are terrible 
fits to a host of astronomical data: large-scale structure observations, supernova data 
and measurements of local dynamics. As discussed in §7.3, the combination of WMAP 
data and other astronomical data solidifies the evidence against these models. The de- 
tected cross-correlation between CMB fluctuations and large-scale structure provides 
further evidence for the existence of dark energy (see §4.1.10). 

• The simple scale invariant (n s = 1.0) model is no longer a good fit to the WMAP 
data. As discussed in the previous subsection, combining the WMAP data with other 
astronomical data sets further strengthens the case for n s < 1. 

The conclusion that the WMAP data demands the existence of dark matter and dark energy 
is based on the assumption that the primordial power spectrum is a power-law spectrum. By 
adding additional features in the primordial perturbation spectrum, these alternative models 
may be able to better mimic the ACDM model. This possibility requires further study. 

The bottom half of Table 3 lists the relative improvement of the generalized models 
over the power-law ACDM. As the Table shows, the WMAP data alone does not require 
the existence of tensor modes, quintessence, or modifications in neutrino properties. Adding 
these parameters does not improve the fit significantly. For the WMAP data, the region in 
likelihood space where (r = 0, w = —1, and m v = lies within the la contour. In the §7, 
we consider the limits on these parameters based on WMAP data and other astronomical 
data sets. 

If we allow for a non-flat universe, then models with small negative Qk are a better 
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fit than the power-law ACDM model. These models have a lower intervening Sachs- Wolfe 
(ISW) signal at low I and are a better fit to the low £ multipoles. The best fit closed 
universe model has Q m = 0.415, = 0.630 and H = 55 kms _1 Mpc _1 and is a better fit 
to the WMAP data alone than the flat universe model(Axlff = 2) However, as discussed 
in §7.3, the combination of WMAP data with either SNe data, large-scale structure data or 
measurements of Hq favors models with Qk close to 0. 

In section 5, we consider several different modifications to the shape of the power spec- 
trum. As noted in Table 3 , none of the modifications lead to significant improvements in the 
fit. Allowing the spectral index to vary as a function of scale improves the goodness-of-fit. 
The low I multipoles, particularly i = 2, are lower than predicted in the ACDM model. 
However, the relative improvement in fit is not large, &xlff = 3, so the WMAP data alone 
do not require a running spectral index. 

Measurement of the goodness of fit is a simple approach to test the needed number of 
parameters. These results should be confirmed by Bayesian model comparison techniques 
(Beltran et al. 2005; Trotta 2005; Mukherjee et al. 2006; Bridges et al. 2005). Bayesian 
methods, however, require an estimate of the number of data points in the fit. It is not clear 
whether we should use the ~ 10 10 points in the TOD, the 10 6 points in the temperature 
maps, the 3 x 10 3 points in the TT, TE, or EE power spectrum, or the ~ 10 — 20 numbers 
needed to fit the peaks and valleys in the TT data in evaluating the significance of new 
parameters. 

4. WMAP ACDM Model and Other Astronomical Data 

In this paper, our approach is to show first that a wide range of astronomical data sets 
are consistent with the predictions of the ACDM model with its parameters fitted to the 
WMAP data (see section §4.1). We then use the external data sets to constrain extensions 
of the standard model. 

In our analyses, we consider several different types of data sets. We consider the SDSS 
LRGs, the SDSS full sample and the 2dFGRS data separately, this allows a check of system- 
atic effects. We divide the small scale CMB data sets into low frequency experiments (CBI, 
VSA) and high frequency experiments (BOOMERanG, ACBAR). We divide the supernova 
data sets into two groups as described below. The details of the data sets are also described 
in §4.1. 

When we consider models with more parameters, there are significant degeneracies, and 
external data sets are essential for parameter constraints. We use this approach in §4.2 and 
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subsequent sections. 

4.1. Predictions from the WMAP Best Fit ACDM Model 

The WMAP data alone is now able to accurately constrain the basic six parameters 
of the ACDM model. In this section, we focus on this model and begin by using only the 
WMAP data to fix the cosmological parameters. We then use the Markov chains (and linear 
theory) to predict the amplitude of fluctuations in the local universe and compare to other 
astronomical observations. These comparisons test the basic physical assumptions of the 
ACDM model. 



4-1.1. Age of the Universe and Hq 

The CMB data do not directly measure H ; however, by measuring il m Ho through the 
height of the peaks and the conformal distance to the surface of last scatter through the peak 
positions (Page et al. 2003b), the CMB data produces a determination of H if we assume the 
simple flat ACDM model. Within the context of the basic model of adiabatic fluctuations, 
the CMB data provides a relatively robust determination of the age of the universe as the 
degeneracy in other cosmological parameters is nearly orthogonal to measurements of the 
age (Knox et al. 2001; Hu et al. 2001). 

The WMAP ACDM best fit value for the age: t = 13.731^ Gyr, agrees with estimates 
of ages based on globular clusters (Chaboyer & Krauss 2002) and white dwarfs (Hansen et al. 
2004; Richer et al. 2004). Figure 4 compares the predicted evolution of H(z) to the HST 
key project value (Freedman et al. 2001) and to values from analysis of differential ages as 
a function of redshift (Jimenez et al. 2003; Simon et al. 2005). The WMAP best fit value, 
H =73.2+3;2 km/s/Mpc, is also consistent with HST measurements (Freedman et al. 2001), 
H = 72 ± 8 km/s/Mpc, where the error includes random and systematic uncertainties and 
the estimate is based on several different methods (Type la supernovae, Type II supernovae, 
surface brightness fluctuations and fundamental plane). It also agrees with detailed studies of 
gravitationally lensed systems such as B 1608+656 (Koopmans et al. 2003), which yields 75TJg 
km/s/Mpc, measurements of the Hubble constant from SZ and X-ray observations of clusters 
(Bonamente et al. 2005) that find H = 76^4 lg° ° km/s/Mpc, and recent measurements 
of the Cepheid distances to nearby galaxies that host type la supernova (Riess et al. 2005), 
H = 73 ± 4 ± 5 km/s/Mpc. 
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Fig. 4. The ACDM model fit to the WMAP data predicts the Hubble parameter 
redshift relation. The blue band shows the 68% confidence interval for the 
Hubble parameter, H. The dark blue rectangle shows the HST key project 
estimate for H and its uncertainties (Freedman et al. 2001). The other points are 
from measurements of the differential ages of galaxies, based on fits of synthetic 
stellar population models to galaxy spectroscopy. The squares show values from 
Jimenez et al. (2003) analyses of SDSS galaxies. The diamonds show values from 
Simon et al. (2005) analysis of a high redshift sample of red galaxies. 
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4-1.2. Big Bang Nucleosynthesis 

Measurements of the light element abundances are among the most important tests of 
the standard big bang model. The WMAP estimate of the baryon abundance depends on our 
understanding of acoustic oscillations 300,000 years after the big bang. The BBN abundance 
predictions depend on our understanding of physics in the first minutes after the big bang. 

Table 4 lists the primordial deuterium abundance, y^ IT , the primordial 3 He abundance, 
ys, the primordial helium abundance, Yp, and the primordial 7 Li abundance, yu, based on 
analytical fits to the predicted BBN abundances (Kneller & Steigman 2004) and the power- 
law ACDM 68% confidence range for the baryon/photon ratio, r] 10 = (273.9 ± 0.3)il b h 2 . The 
lithium abundance is often expressed as a logarithmic abundance, [Li]p = 12 + logio(Li/H). 

Table 4: Primordial abundances based on using the Steigman (2005) fitting formula for the 
ACDM 3 -year WMAP only value for the baryon/photon ratio, t]iq — 6.116^q^49- 





CMB-based BBN prediction 


Observed Value 




9 p-7+0.17 
' -0.13 


1.6 - 4.0 


10 5 t/3 


1.05 ±0.03 ±0.03 (syst.) 


< 1.1 ±0.2 


Yp 


0.24819lg;K ± 0.0006(syst.) 


0.232 - 0.258 


[Li] 


2.64 ±0.03 


2.2 - 2.4 



The systematic uncertainties in the helium abundances are due to the uncertainties in 
nuclear parameters, particularly neutron lifetime (Steigman 2005). Prior to the measure- 
ments of the CMB power spectrum, uncertainties in the baryon abundance were the biggest 
source of uncertainty in CMB predictions. Serebrov et al. (2005) argues that the currently 
accepted value, r„ = 887.5 s, should be reduced by 7.2 s, a shift of several times the reported 
errors in the Particle Data Book. This (controversial) shorter lifetime lowers the predicted 
best fit helium abundance to Y P = 0.24675 (Mathews et al. 2005; Steigman 2005). 

The deuterium abundance measurements provide the strongest test of the predicted 
baryon abundance. Kirkman et al. (2003) estimate a primordial deuterium abundance, 
[D]/[H]= 2.78+o; 38 x 10~ 5 , based on five QSO absorption systems. The six systems used in 
the Kirkman et al. (2003) analysis show a significant range in abundances: 1.65 — 3.98 x 10~ 5 
and have a scatter much larger than the quoted observational errors. Recently, Crighton et al. 
(2004) report a deuterium abundance of 1.6l°;4 x 10 ~ 5 for PKS 1937-1009. Because of the 
large scatter, we quote the range in [D]/[H] abundances in Table 4; however, note that the 
mean abundance is in good agreement with the CMB prediction. 
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It is difficult to directly measure the primordial 3 He abundance. Bania et al. (2002) 
quote an upper limit on the primordial 3 He abundance of 2/3 < 1.1 ± 0.2 x 10~ 5 . This limit 
is compatible with the BBN predictions. 

Olive & Skillman (2004) have reanalyzed the estimates of primordial helium abundance 
based on observations of metal-poor HII regions. They conclude that the errors in the 
abundance are dominated by systematic errors and argue that a number of these systematic 
effects have not been properly included in earlier analysis. In Table 4, we quote their estimate 
of the allowed range of Yp. Olive & Skillman (2004) find a representative value of 0.249±0.009 
for a linear fit of [0]/[H] to the helium abundance, significantly above earlier estimates and 
consistent with WMAP- normalized BBN. 

While the measured abundances of the other light elements appear to be consistent 
with BBN predictions, measurements of neutral lithium abundance in low metallicity stars 
imply values that are a factor of 2 below the BBN predictions: most recent measurements 
(Charbonnel & Primas 2005; Boesgaard et al. 2005) find an abundance of [Li]p ~ 2.2 — 2.25. 
While Melendez & Ramirez (2004) find a higher value, [Li] P ~ 2.37 ± 0.05, even this value 
is still significantly below the BB-predicted value, 2.64 ±0.03. These discrepancies could be 
due to systematics in the inferred lithium abundance (Steigman 2005), uncertainties in the 
stellar temperature scale (Fields et al. 2005), destruction of lithium in an early generation 
of stars or the signature of new early universe physics (Coc et al. 2004; Jedamzik 2004; 
Richard et al. 2005; Ellis et al. 2005; Larena et al. 2005; Jedamzik et al. 2005). The recent 
detection (Asplund et al. 2005) of 6 Li in several low metallicity stars further constrains 
chemical evolution models and exacerbates the tensions between the BBN predictions and 
observations. 



4-1.3. Small Scale CMB Measurements 

With the parameters of the A CDM model fixed by the measurements of the first three 
acoustic peaks, the basic properties of the small scale CMB fluctuations are determined by the 
assumption of a power law for the amplitude of potential fluctuations and by the physics of 
Silk damping. We test these assumptions by comparing the WMAP best fit power law ACDM 
model to data from several recent small scale CMB experiments (BOOMERanG, MAXIMA, 
ACBAR, CBI, VSA). These experiments probe smaller angular scales than the WMAP 
experiment and are more sensitive to the details of recombination and the physics of acoustic 
oscillations. The good agreement seen in Figure (5) suggests that the standard cosmological 
model is accurately characterizing the basic physics at z ~ 1100. In subsequent sections, we 
combine WMAP with small scale experiments. We include four external CMB datasets which 
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Fig. 5. — The prediction for the small-scale angular power spectrum seen by 
ground-based and balloon CMB experiments from the ACDM model fit to the 
WMAP data only. The colored lines show the best fit (red) and the 68% (dark 
orange) and 95% confidence levels (light orange) based on fits of the ACDM 
models to the WMAP data. The points in the figure show small scale CMB 
measurements (Ruhl et al. 2003; Abroe et al. 2004; Kuo et al. 2004; Readhead 
et al. 2004a; Dickinson et al. 2004). The plot shows that the ACDM model (fit 
to the WMAP data alone) can accurately predict the amplitude of fluctuations 
on the small scales measured by ground and balloon-based experiments. 

complement the WMAP observations at smaller scales: the Cosmic Background Imager 
(CBI: Mason et al. (2003); Sievers et al. (2003); Pearson et al. (2003); Readhead et al. 
(2004a)), the Very Small Array (VSA: Dickinson et al. (2004)), the Arcminute Cosmology 
Bolometer Array Receiver (ACBAR: Kuo et al. (2004)) and BOOMERanG (Ruhl et al. 
2003; Montroy et al. 2005; Piacentini et al. 2005) We do not include results from a number 
of experiments that overlap in i range coverage with WMAP as these experiments have 
non-trivial cross-correlations with WMAP that would have to be included in the analysis. 

We do not use the small-scale polarization results for parameter determination as they do 
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not yet noticeably improve constraints. These polarization measurements, however, already 
provide important tests on the basic assumptions of the model (e.g., adiabatic fluctuations 
and standard recombination history). 

The measurements beyond the third peak improve constraints on the cosmological pa- 
rameters. These observations constrict the {r, Ub, A s ,n s } degeneracy and provide an im- 
proved probe of a running tilt in the primordial power spectrum. In each case we only use 
bandpowers that do not overlap with signal-dominated WMAP observations, so that they 
can be treated as independent measurements. 

In the subsequent sections, we perform likelihood analysis for two combinations of 
WMAP data with other CMB data sets: WMAP + high frequency bolometer experi- 
ments (ACBAR + BOOMERanG) and WMAP + low frequency radiometer experiments 
(CBI+VSA). The CBI data set is described in Readhead et al. (2004a). We use 7 band pow- 
ers, with mean I values of 948, 1066, 1211, 1355, 1482, 1692 and 1739, from the even binning 
of observations rescaled to a 32 GHz Jupiter temperature of 147.3 ± 1.8K. The rescaling 
reduces the calibration uncertainty to 2.6% from 10% assumed in the first year analyses; 
CBI beam uncertainties scale the entire power spectrum and, thus, act like a calibration 
error. We use a log- normal form of the likelihood as in Pearson et al. (2003). The VSA 
data (Dickinson et al. 2004) uses 5 band powers with mean t- values of 894, 995, 1117, 1269 
and 1407, which are calibrated to the WMAP 32 GHz Jupiter temperature measurement. 
The calibration uncertainty is assumed to be 3% and again we use a log-normal form of the 
likelihood. For ACBAR (Kuo et al. 2004), we use the same band powers as in the first year 
analysis, with central t values 842, 986, 1128, 1279, 1426, 1580, and 1716, and errors from 
the ACBAR web site 1 . We assume a calibration uncertainty of 20% in Ce, and the quoted 
3% beam uncertainty in Full Width Half Maximum. We use the temperature data from the 
2003 flight of BOOMERanG, based on the "NA pipeline" (Jones et al. 2005) considering 
the 7 datapoints and covariance matrix for bins with mean £ values, 924, 974, 1025, 1076, 
1117, 1211 and 1370. While there is overlap between some of the small scale power spectrum 
measurements and WMAP measurements at 800 < I < 1000, the WMAP measurements are 
noise-dominated so that there is little covariance between the measurements. 



4-1. 4- Large-Scale Structure 

With the WMAP polarization measurements constraining the suppression of temper- 
ature anisotropy by reionization, we now have an accurate measure of the amplitude of 



See http://cosmology.berkeley.edu/group/swlh/acbar/data 
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fluctuations at the redshift of decoupling. If the power-law ACDM model is an accurate 
description of the large-scale properties of the universe, then we can extrapolate forward 
the roughly 1000-fold growth in the amplitude of fluctuations due to gravitational clustering 
and compare the predicted amplitude of fluctuations to the large-scale structure observa- 
tions. This is a powerful test of the theory as some alternative models fit the CMB data but 
predict significantly different galaxy power spectra (e.g., Blanchard et al. (2003)). 

Using only the WMAP data, together with linear theory, we can predict the amplitude 
and shape of the matter power spectrum. The band in Figure 6 shows the 68% confidence 
interval for the matter power spectrum. Most of the uncertainty in the figure is due to 
the uncertainties in Q m h. The points in the figure show the SDSS Galaxy power spectrum 
(Tegmark et al. 2004b) with the amplitude of the fluctuations normalized by the galaxy 
lensing measurements and the 2dFGRS data (Cole et al. 2005). The figure shows that the 
ACDM model, when normalized to observations at z ~ 1100, accurately predicts the large- 
scale properties of the matter distribution in the nearby universe. It also shows that adding 
the large-scale structure measurements will reduce uncertainties in cosmological parameters. 




Fig. 6. — (Left) The predicted power spectrum (based on the range of parameters 
consistent with the WMAP-only parameters) is compared to the mass power 
spectrum inferred from the SDSS galaxy power spectrum (Tegmark et al. 2004b) 
as normalized by weak lensing measurements (Seljak et al. 2005b). (Right) The 
predicted power spectrum is compared to the mass power spectrum inferred 
from the 2dFGRS galaxy power spectrum (Cole et al. 2005) with the best fit 
value for b 2d FGRS based on the fit to the WMAP model. Note that the 2dFGRS 
data points shown are correlated. 

When we combine WMAP with large-scale structure observations in subsequent sections, 
we consider the measurements of the power spectrum from the two large-scale structure sur- 
veys. Since the galaxy power spectrum does not suffer the optical depth-driven suppression 
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in power seen in the CMB, large scale structure data gives an independent measure of the 
normalization of the initial power spectrum (to within the uncertainty of the galaxy biasing 
and redshift space distortions) and significantly truncates the {r, u;&, A s ,n s } degeneracy. In 
addition the galaxy power spectrum shape is determined by Q m h as opposed to the fl m h 2 
dependency of the CMB. Its inclusion therefore further helps to break the {cu m , w or Qk} 
degeneracy. 

The 2dFGRS survey probes the universe at redshift z ~0.1 (we assume z e ff = 0.17 for 
the effective redshift for the survey) and probes the power spectrum on scales 0.022 /iMpc -1 < 
k < 0.19 /iMpc -1 . From the data and covariance described in Cole et al. (2005) we use 32 
of the 36 bandpowers in the range 0.022 /iMpc -1 < k < 0.19 /iMpc -1 . We correct for non- 
linearities and non-linear redshift space distortions using the prescription employed by the 
2dF team, 

p r 9 :t sh {k) = 1 - T ^p t ; a r(k) (5) 

where P^j* 1 and Pg^f° ry are the redshift space and theoretical real space galaxy power 
spectra with Q = 4. (Mpc/h) 2 and A = 1.4 Mpc/h. We analytically marginalize over the 
power spectrum amplitude, effectively applying no prior on the linear bias and on linear 
redshift space distortions, in contrast to our first year analyses. 

The SDSS main galaxy survey measures the galaxy distribution at redshift of z ~ 0.1; 
however, as in the analysis of the SDSS team (Tegmark et al. 2004b) we assume z e ff = 0.07, 
and we use 14 power spectrum bandpowers between O.OI6/1 Mpc -1 < k < 0.11/t Mpc -1 . 
We follow the approach used in the SDSS analysis in Tegmark et al. (2004a): we use the 
nonlinear evolution of clustering as described in Smith et al. (2003) and include a linear bias 
factor, b s dss, and the linear redshift space distortion parameter, (5. 

PiTik) = (1 + \{3 + \f)Pf a r V {k) (6) 

Following Lahav et al. (1991), we use (3b = d In 5/d In a where (3 w [fi^ 7 +(l+fi m /2)(fi A /70)]/6. 
We impose a Gaussian prior on the bias of b = 1.03 +/- 0.15, based on an estimate from 
weak lensing of the same SDSS galaxies used to derive the matter power spectrum. This 
value includes a 4% calibration uncertainty in quadrature with the reported bias error 2 and 
is a symmetrized form of the bias constraint in Table 2 of Seljak et al. (2005b). While the 
WMAP first year data was used in the Seljak et al. (2005b) analysis, the covariance between 
the data sets are small. We restrict our analysis to scales where the bias of a given galaxy 
population does not show significant scale dependence (Zehavi et al. 2005). Analyses that 



2 M. Tegmark private communication. 
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use galaxy clustering data on smaller scales require detailed modeling of non-linear dynamics 
and the relationship between galaxy halos and galaxy properties (see, e.g., Abazajian et al. 
(2005)). 

The SDSS luminous red galaxy (LRG) survey uses the brightest class of galaxies in the 
SDSS survey (Eisenstein et al. 2005). While a much smaller galaxy sample than the main 
SDSS galaxy sample, it has the advantage of covering 0.72 h~ 3 Gpc 3 between 0.16 < z < 0.47. 
Because of its large volume, this survey was able to detect the acoustic peak in the galaxy 
correlation, one of the distinctive predictions of the standard adiabatic cosmological model 
(Peebles & Yu 1970; Sunyaev & Zel'dovich 1970; Silk & Wilson 1981; Bond & Efstathiou 
1984; Vittorio & Silk 1984; Bond & Efstathiou 1987). We use the SDSS acoustic peak 
results to constrain the balance of the matter content, using the well measured combination, 



A(z = 0.35) 



d A (z = 0.35) 2 
H(z = 0.35) 



(7) 



where oIa is the comoving angular diameter distance and c is the speed of light. We impose 
a Gaussian prior of A = 0.469 (^fg) ° 35 ± 0.017 based on the analysis of Eisenstein et al. 
(2005) . 



4-1.5. Lyman a Forest 

Absorption features in high redshift quasars (QSO) at around the frequency of the 
Lyman-a emission line are thought to be produced by regions of low-density gas at redshifts 
2 < z < 4 (Croft et al. 1998; Gnedin & Hamilton 2002). These features allow the matter 
distribution to be characterized on scales of 0.2 < k < 5 h Mpc -1 and as such extend the 
lever arm provided by combining large-scale structure data and CMB. These observations 
also probe a higher redshift range (z ~ 2 — 3). Thus, these observations nicely complement 
CMB measurements and large scale structure observations. While there has been significant 
progress in understanding systematics in the past few years (McDonald et al. 2005; Meiksin 
& White 2004), time constraints limit our ability to consider all relevant data sets. 

Recent fits to the Lyman-a forest imply a higher amplitude of density fluctuations: Jena 
et al. (2005) find that <r 8 = 0.9, tt m = 0.27, h = 0.71 provides a good fit to the Lyman a data. 
Seljak et al. (2005a) combines first year WMAP data, other CMB experiments, large-scale 
structure and Lyman a to find: n s = 0.98 ± 0.02, a 8 = 0.90 ± 0.03, h = 0.71 ± 0.021, and 
Q m = 0.281+o;q2?- Note that if they assume r = 0.09, the best fit value drops to a 8 = 0.84. 
While these models have somewhat higher amplitudes than the new best fit WMAP values, 
a recent analysis by Desjacques & Nusser (2005) find that the Lyman a data is consistent 
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with as between 0.7 — 0.9. This suggests that the Lyman a data is consistent with the new 
WMAP best fit values; however, further analysis is needed. 

4-1.6. Galaxy Motions and Properties 

Observations of galaxy peculiar velocities probe the growth rate of structure and are 
sensitive to the matter density and the amplitude of mass fluctuations. The Feldman et al. 
(2003) analysis of peculiar velocities of nearby ellipticals and spirals finds Vt m = 0.30±o;o7 
and as = 1.13lo!23> within la of the WMAP best fit value for Vt m and 1.5a higher than the 
WMAP value for as- These estimates are based on dynamics and are not sensitive to the 
shape of the power spectrum. Mohayaee & Tully (2005) apply orbit retracing methods to 
motions in the local supercluster and obtain Q m = 0.22 ± 0.02, consistent with the WMAP 
values. 

Modeled galaxy properties can be compared to the clustering properties of galaxies 
on smaller scales. The best fit parameters for WMAP only are consistent with the recent 
Abazajian et al. (2005) analysis of the pre-three year release CMB data combined with 
the SDSS data. In their analysis, they fit a Halo Occupation Distribution model to the 
galaxy distribution so as to use the galaxy clustering data at smaller scales. Their best fit 
parameters (H = 70 ± 2.6 km/s/Mpc, f2 m = 0.271 ± 0.026) are consistent with the results 
found here. Vale & Ostriker (2005) fit the observed galaxy luminosity functions with a s = 0.8 
and Q m = 0.25. Van den Bosch, Mo & Yang (2003) used the conditional luminosity function 
to fit the 2dFGRS luminosity function and the correlation length as a function of luminosity. 
Combining with the first-year WMAP data, they found Q m = 0.2518;$ and a 8 = 0.78 ±0.12 
(95% CL), again in remarkable agreement with the three year WMAP best fit values. 

4-1.7. Weak Lensing 

Over the past few years, there has been dramatic progress in using weak lensing data 
as a probe of mass fluctuations in the nearby universe. Lensing surveys complement CMB 
measurements (Contaldi et al. 2003; Tereno et al. 2005), and their dominant systematic 
uncertainties differ from the large-scale structure surveys. 

Measurements of weak gravitational lensing, the distortion of galaxy images by the 
distribution of mass along the line of sight, directly probe the distribution of mass fluctuations 
along the line of sight (see Refregier (2003) for a recent review). Figure 7 shows that the 
WMAP ACDM model predictions for as and Q m are lower than the amplitude found in 
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Fig. 7. — Prediction for the mass fluctuations measured by the CFTHLS weak- 
lensing survey from the ACDM model fit to the WMAP data only. The blue, 
red and green contours show the joint 2-d marginalized 68% and 95% confidence 
limits in the (a s , tt m ) plane for for WMAP only, CFHTLS only and WMAP 
+ CFHTLS, respectively, for the power law ACDM models. All constraints 
come from assuming the same priors on input parameters, with the additional 
marginalization over z s in the weak lensing analysis, using a top hat prior of 
0.613 < z s < 0.721 . While lensing data favors higher values of a 8 ~ 0.8 — 1.0 (see 
§4.1.7), X-ray cluster studies favor lower values of as — 0.7 — 0.8 (see §4.1.9). 
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most recent lensing surveys: Hoekstra et al. (2002) calculate a% = 0.94+Q;}4(f2 m /0.25) °- 52 
(95% confidence) from the RCS survey and Van Waerbeke et al. (2005) determine as = 
0.91±0.08(fi m /0.25)-° 49 from the VIRMOS-DESCART survey; however, Jarvis et al. (2003) 
find <t 8 = 0.79^;^(fi m /0.25)-°- 57 (95% confidence level) from the 75 Degree CTIO survey. 

In §4.2, we use the data set provided by the first weak gravitational lensing analysis of 
the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) 3 as conducted by Hoekstra 
et al. (2005) (Ho05) and Semboloni et al. (2005). Following Ho05, we use only the wide fields 
Wl and W3, hence a total area of 22 deg 2 observed in the i' band limited to a magnitude of 
i' = 24.5. We follow the same methodology as Ho05 and Tereno et al. (2005). For each given 
model and set of parameters, we compute the predicted shear variance at various smoothing 
scales, (7 2 ), and then evaluate its likelihood (see Ho05 equation 13). 

Since the lensing data are in a noise dominated regime, we neglect the cosmological 
dependence of the covariance matrix. To account conservatively for a possible residual 
systematic contamination, we use (7^) as a monitor and add it in quadrature to the diagonal 
of the noise covariance matrix, as in Ho05. We furthermore marginalize over the mean source 
redshift, z s (defined in equation (16) of Ho05) assuming a uniform prior between 0.613 and 
0.721. This marginalization is performed by including these extra parameters in the Monte 
Carlo Markov Chain. Our analysis differs however from the likelihood analysis of Ho05 
in the choice of the transfer function. We use the Novosyadlyj et al. (1999) (NDL) CDM 
transfer function (with the assumptions of Tegmark et al. (2001)) rather than the Bardeen 
et al. (1986) (BBKS) CDM transfer function. The NDL transfer function includes more 
accurately baryon oscillations and neutrino effects. This modification alters the shape of the 
likelihood surface in the 2-dimensional (<7 8 ,f2 m ) likelihood space. 



4-1.8. Strong Lensing 

Strong lensing provides another potentially powerful probe of cosmology. The number 
of multiply-lensed arcs and quasars is very sensitive to the underlying cosmology. The 
cross-section for lensing depends on the number of systems with surface densities above 
the critical density, which in turn is sensitive to the angular diameter distance relation 
(Turner 1990). The CLASS lensing survey (Chae et al. 2002) finds that the number of lenses 
detected in the radio survey is consistent with a flat universe with a cosmological constant 
and tt m = 0.31+g;fl. The statistics of strong lenses in the SDSS is also consistent with the 
standard ACDM cosmology (Oguri 2004). The number and the properties of lensed arcs are 



3 http://www.cfht.hawaii.edu/Science/CFHTLS 
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also quite sensitive to cosmological parameters (but also to the details of the data analysis). 
Wambsganss et al. (2004) conclude that arc statistics are consistent with the concordance 
ACDM model. 

Soucail et al. (2004) has used multiple lenses in Abell 2218 to provide another geomet- 
rical tests of cosmological parameters. They find that < il m < 0.33 and w < —0.85 for 
a flat universe with dark energy. This method is another independent test of the standard 
cosmology. 

4-1.9. Clusters and the Growth of Structure 

The numbers and properties of rich clusters are another tool for testing the emerging 
standard model. Since clusters are rare, the number of clusters as a function of redshift is 
a sensitive probe of cosmological parameters. Recent analyses of both optical and X-ray 
cluster samples yield cosmological parameters consistent with the best fit WMAP ACDM 
model (Borgani et al. 2001; Bahcall & Bode 2003; Allen et al. 2003; Vikhlinin et al. 2003; 
Henry 2004). The parameters are, however, sensitive to uncertainties in the conversion 
between observed properties and cluster mass (Pierpaoli et al. 2003; Rasia et al. 2005). 

Clusters can also be used to infer cosmological parameters through measurements of the 
baryon/dark matter ratio as a function of redshift (Pen 1997; Ettori et al. 2003; Allen et al. 
2004). Under the assumption that the baryon/dark matter ratio is constant with redshift, 
the Universe is flat, and standard baryon densities, Allen et al. (2004) find Vt m = 0.24 ±0.04 
and w = -1.20^. Voevodkin & Vikhlinin (2004) determine cr 8 = 0.72 ± 0.04 and Q m h 2 = 
0.13 ±0.07 from measurements of the baryon fraction. These parameters are consistent with 
the values found here and in §7.1. 

4.1.10. Integrated Sachs- Wolfe (IS W) effect 

The ACDM model predicts a statistical correlation between CMB temperature fluctua- 
tions and the large-scale distribution of matter (Crittenden & Turok 1996). Several groups 
have detected correlations between the WMAP measurements and various tracers of large- 
scale structure at levels consistent with the concordance ACDM model (Boughn & Crittenden 
2004; Nolta et al. 2004; Afshordi et al. 2004; Scranton et al. 2003; Fosalba k Gaztahaga 2004; 
Padmanabhan et al. 2005; Corasaniti et al. 2005; Boughn & Crittenden 2005; Vielva et al. 
2006). These detections provide an important independent test of the effects of dark energy 
on the growth of structure. However, the first year WMAP data is already signal dominated 
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on the scales probed by the ISW effect, thus, improved large-scale structure surveys are 
needed to improve the statistical significance of this detection (Afshordi 2004; Bean & Dore 
2004; Pogosian et al. 2005). 



4-1.11. Supernova 




Fig. 8. — Using the WMAP A CDM parameters, we predict luminosity distance- 
redshift relationship and compare it to measurements from supernova surveys. 
The plots show the deviations of the distance measure (DM) from the empty 
universe model. The solid lines are for the best WMAP ACDM parameters and 
the blue band shows the 68% confidence range. (Left) SNLS data(Astier et al. 
2005). (Right) "gold" supernova data (Riess et al. 2004). 

With the realization that their light curve shapes could be used to make SN la into 
standard candles, supernovae have become an important cosmological probe (Phillips 1993; 
Hamuy et al. 1996; Riess et al. 1996). They can be used to measure the luminosity distance 
as a function of redshift. The dimness of z ~ 0.5 supernova provide direct evidence for 
the accelerating universe (Riess et al. 1998; Schmidt et al. 1998; Perlmutter et al. 1999; 
Tonry et al. 2003; Knop et al. 2003; Nobili et al. 2005; Clocchiatti et al. 2005; Krisciunas 
et al. 2005; Astier et al. 2005). Recent HST measurements (Riess et al. 2004) trace the 
luminosity distance/redshift relation out to higher redshift and provide additional evidence 
for presence of dark energy. Assuming a flat Universe, the Riess et al. (2004) analysis of 
the supernova data alone finds that Vt m = 0.29loo3 consistent with the fits to WMAP data 
alone (see Table 2) and to various combinations of CMB and LSS data sets (see Tables 5 
and 6). Astier et al. (2005) find that Q m = 0.263t° o fi(stat.)toml(sys.) from the first year 
supernova legacy survey. 

Within the ACDM model, the supernovae data serve as a test of our cosmological model. 



-33- 



Figure 8 shows the consistency between the supernova and CMB data. Using just the WMAP 
data and the ACDM model, we can predict the distance/luminosity relationship and test it 
with the supernova data. 

In §4.2 and subsequent sections, we consider two recently published high-z supernovae 
datasets in combination with the WMAP CMB data, the first sample is 157 supernova 
in the "Gold Sample" as described in Riess et al. (2004) with 0.015 < z < 1.6 based on a 
combination of ground-based data and the GOODS ACS Treasury program using the Hubble 
Space Telescope (HST) and the second sample is 115 supernova in the range 0.015 < z < 1 
from the Supernova Legacy Survey (SNLS) (Astier et al. 2005) . 

Measurements of the apparent magnitude, m, and inferred absolute magnitude, M , 
of each SN have been used to derive the distance modulus fi b s = m — M , from which 
a luminosity distance is inferred, fi b s = 5 log[d l(z) /Mpc] + 25. The luminosity distance 
predicted from theory, fj,th, is compared to observations using a x 2 analysis summing over 
the SN sample. 



where the absolute magnitude, M , is a "nuisance parameter", analytically marginalized over 
in the likelihood analysis (Lewis & Bridle 2002), and a obs contains systematic errors related 
to the light curve stretch factor, K-correction, extinction and the intrinsic redshift dispersion 
due to SNe peculiar velocities (assumed 400 and 300 km s" 1 for HST/GOODS and SNLS 
data sets respectively). 



In the previous section we showed that the power-law LCDM model fit to WMAP 
data only is consistent with other astronomical data. Motivated by this stringent series of 
cosmological tests, we combine the WMAP data with other astronomical observations to 
further constrain cosmological parameters. 

Tables 5 and 6 show that adding external data sets has little effect on several parameters: 
Qbh 2 , n s and r. However, the various combinations do reduce the uncertainties on Q m and 
the amplitude of fluctuations. The data sets used in Table 5 favor smaller values of the 
matter density, higher Hubble constant values, and lower values of a 8 . The data sets used 
in Table 6 favor higher values of Q m , lower Hubble constants and higher values of a 8 . The 
lensing data is most discrepant and it most strongly pulls the combined results towards 
higher amplitudes and higher Q m (see Figure 7 and 9). The overall effect of combining the 




(8) 



4.2. 



Joint Constraints on ACDM Model Parameters 
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Table 5: ACDM Model: Joint Likelihoods These values are calculated using the N side = 8 likelihood 
code with Apg = 0.017 
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Table 6: ACDM Model 
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data sets is shown in Figure 10. 

The best fits shown in Table 6 differ by about la from the best fits shown in Table 5 
in their predictions for the total matter density, Q m h 2 (See Tables 5 and 6 and Figure 9). 
More accurate measurements of the third peak will help resolve these discrepancies. 

The differences between the two sets of data may be due to statistical fluctuations. 
For example, the SDSS main galaxy sample power spectrum differs from the power spec- 
trum measured from the 2dfGRS: this leads to a lower value for the Hubble constant for 
WMAP±SDSS data combination, h = 0.710 ± 0.026, than for WMAP±2dFGRS, h = 
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0.733_ 021 . Note that while the SDSS LRG data parameters values are close to those 
from the main SDSS catalog, they are independent determinations with mostly different 
systematics. 

Lensing measurements are sensitive to the amplitude of the local potential fluctuations, 
which scale roughly as <7gQ^, so that lensing parameter constraints are nearly orthogonal 
to the CMB degeneracies (Tereno et al. 2005). The CFHTLS lensing data best fit value for 
o"gfi ? n 6 is 1 — 2cr higher than the best fit three year WMAP value. As a result, the combination 
of CFHT and WMAP data favors a higher value of erg and Q m and a lower value of H Q than 
WMAP data alone. Appendix A shows that the amplitude of this discrepancy is somewhat 
sensitive to our choice of priors. Because of the small error bars in the CFHT data set 
and the relatively small overlap region in parameter space, the CFHT data set has a strong 
influence on cosmological parameters. Because of the small errors in the CFHT data and 
the relatively small overlap region in parameter space, the CFHT data has a strong influence 
on cosmological parameters. This effect if exacerbated when additional cosmological data 
sets are included in the analysis. Because of this, we do not include the lensing data in the 
full combined data set (WMAP... WMAP+ ALL), rather we quote WMAP+CFHT results 




0.10 0.11 0.12 0.13 0.14 0.15 

m /i 2 

Fig. 9. One-dimensional marginalized distribution of Vt m h 2 for 

WMAP, WMAP+CBI+VSA, WMAP+BOOM+ACBAR, WMAP+SDSS, 
WMAP+SN(SNLS), WMAP+SN(HST/GOODS), WMAP+2dFGRS and 
WMAP+CFHTLS for the power-law ACDM model. 
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Fig. 10. — Joint two-dimensional marginalized contours (68%, and 95% con- 
fidence levels) for various combination of parameters for WMAP only (solid 
lines) and WMAP+2dFGRS+SDSS+ACBAR+BOOMERanG+CBI+VSA+ 
SN(HST/GOODS)+SN(SNLS) (filled red) for the power-law ACDM model. 
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separately. The combined data sets place the strongest limits on cosmological parameters. 
Because they are based on the overlap between many likelihood functions, limits based on 
the WMAP+ALL data set should be treated with some caution. Figure 10 shows the 2- 
dimensional marginalized likelihood surface for both WMAP only and for the combination 
of WMAP+ALL. 



While the simplest inflationary models predict that the spectral index varies slowly 
with scale, inflationary models can produce strong scale dependent fluctuations (see e.g., 
Kawasaki, Yamaguchi & Yokoyama (2003); Hall et al. (2004); Yamaguchi & Yokoyama 
(2004)). The first year WMAP observations provided some motivation for considering these 
models as the data, particularly when combined with the Lyman a forest measurements, 
were better fit by models with a running spectral index (Spergel et al. 2003). Small scale 
CMB measurements (Readhead et al. 2004a) also favor running spectral index models over 
power law models. 

Here, we consider whether a more general function for the primordial power spectrum 
could better fit the new WMAP data. We consider three forms for the power spectrum: 

• A^( k) with a running spectral index: 1+cHn A 2 R (k)/d\n k = n(k )+dn s /d\n(k)\n(k/k ) 

• A^(k) allowed to freely vary in 15 bins in /c-space, with k\ = 0, k 2 = 0.001/Mpc, k 15 = 
0.15/Mpc, k i+ i = lAlki for 3 < i < 14. A^(/c) is given by linear interpolation within 
the bins and A^(k) = A^(0.15/Mpc) for k > 0.15/Mpc. 

• A^(/c) with a sharp k cut off at k = k c , 



Figure 11 shows how WMAP data alone can be used to reconstruct the primordial 
power spectrum as a function of scale, parameterized by logarithmically spaced bins out to 
k = 0.15 Mpc -1 . Even allowing for these additional degrees of freedom, the data prefer a 
nearly featureless power-law power spectrum. 



5. Constraining the Shape of the Primordial Power Spectrum 



A 2 




(9) 



oc 



The deviation of the primordial power spectrum from a simple power law can be most 
simply characterized by a sharp cut-off in the primordial spectrum. Analysis of this model 
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Fig. 11. — The reconstructed primordial curvature fluctuation spectrum, Aj^(k), 
for a ACDM cosmology, in logarithmically spaced k bins, where k is in Mpc" 1 . 
The errors show the 68% (red) and 95% (orange) constraints and the black 
diamonds the peak likelihood value. The dashed line show the values for k = 0. 
Consistent with the predictions of simple inflationary theories, there are no 
significant features in the spectrum. The data are consistent with a nearly scale- 
invariant spectrum. 

finds that putting in a cut off of k c ~ 3 x 10~ 4 /Mpc improves the fit by A% 2 = 1.2, not 
enough to justify a radical change in the primordial spectrum. 

Table 3 demonstrates that, while models with reduced large scale power provide slightly 
improved fits to the CMB data, the improvements do not warrant additional parameters. 



The uncertainties in the measurements of running is slightly improved by including 
the small scale experiments. For models with only scalar fluctuations, the marginalized 
value for the derivative of the spectral index is dn s /d\nk = — 0.055^q;q3° for WMAP only, 



dn s /d\nk = -0.066^;°^ for the WMAP+CBI+VSA data and dnjdlnk = -0.058 ± 0.029 
for WMAP+BOOM+ACBAR. For models with tensors, dn s /d\nk = -0.085 ± 0.043 for 
WMAP only, dn s /d\nk = -0.090^1 for WMAP+CBI+VSA, and dn s /d\nk = -0.082 ± 
0.040 for WMAP+BOOM+ACBAR. As Figure 12 shows, models with negative running 
allow large tensor amplitudes; thus, when we marginalize over r with a flat prior, these 
models favor a more negative running. 



5.1. External Data Sets and the Running Spectral Index 
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Fig. 12. — Joint two-dimensional marginalized contours (68% and 95%) for infla- 
tionary parameters, (r,n s ) (left panel) and (r,dn s /d\nk) (right panel), for Model 
Mil in Table 3, with parameters defined at k = 0.002 Mpc -1 . (Upper) WMAP 
only. (Middle) WMAP+SDSS. (Bottom) WMAP+CBI+VSA. Note that n s > 1 
is favored because r and n s are defined at k — 0.002 Mpc -1 . At k = 0.05 Mpc -1 
n s < 1 is favored. The data do not require a running spectral index, dn s /dlnk, at 
more than the 95% confidence level. 



Figure 13 shows that both the power law ACDM model and the running spectral index 
model fit the CMB data. At present, the small scale data do not yet clearly distinguish the 
two models. 

A large absolute value of running would be problematic for most inflationary mod- 
els, so further testing of this suggestive trend is important for our understanding of early 
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Fig. 13. — The running spectral index model provides a slightly better fit to 
the data than the power-law spectral index model. The solid line shows the 
best fit power law ACDM model and the dashed line shows the best fit running 
spectral index ACDM model (fit to WMAP+CBI+VSA). The insert compares 
the models to the WMAP t < 20 data and shows that the running spectral index 
model better fits the decline at £ = 2; however, the improvement in \ 2 ls only 3, 
not enough to strongly argue for the addition of a new parameter. We have also 
done the same analysis for BOOMERanG and ACBAR data and found similar 
results: the current high I data are not yet able to distinguish between the 
running spectral index and power law models. 

universe physics. Additional WMAP data and upcoming, more sensitive small-scale CMB 
data will further test deviations from scale invariance. Current large scale structure data 
do not strengthen the case for running because these these data sets probe similar physical 
scales to the WMAP experiment. Figure 12 shows that current data favor a large running; 
however, the evidence is not yet compelling. The constraints from WMAP+lensing and 
WMAP+2dFGRS are similar to the WMAP+SDSS constraints shown in Figure 12. 



-41 - 



5.2. Is the Power Spectrum Featureless? 

Since inflation magnifies fluctuations that were once on sub-Planckian scales to scales 
of the observable horizon, trans-Planckian physics could potentially leave its imprint on the 
CMB sky. Over the past few years, there has been significant interest in the possibility of 
detecting the signature of trans-Planckian physics in the angular power spectrum. Several 
studies (Martin & Brandenberger 2001; Danielsson 2002; Easther et al. 2002; Bergstrom & 
Danielsson 2002; Kaloper et al. 2002; Martin & Brandenberger 2003; Martin & Ringeval 
2004; Burgess et al. 2003; Schalm, Shiu & van den Schaar 2004) have discussed the possible 
form and the expected amplitude of the trans-Planckian effects which might modify the 
spectrum coming from slow roll inflation. The scalar power spectra resulting from power law 
(PL) slow roll inflation can be written in the terms of Hubble Flow parameters respectively 
as, (Leach & Liddle 2003) 

A 2 n PL (k) =A s (l-2(C+ l)ei - Ce 2 - (2 Cl + e 2 ) In ( f) ) (10) 



k, 



Here, ei and e 2 are slow roll parameters (Leach & Liddle 2003). After the release of the first 
year WMAP data, Martin & Ringeval (2004) considered a primordial power spectrum of a 
slightly modified form to account for additional trans-Planckian (TP) features, 



A^ TP (A;) = A^ PL (k) [1 — 2\x\a cos6>(A;)] — A s \x\a ir(2e 1 + e 2 ) sin9(k) 



with, 9(k) = ^ (l + e 1 + ei In ( 



;n) 



Here a = Hl c /2ir is determined by the Hubble parameter during inflation, H, and the 
characteristic length scale for the trans-Planckian manifestation l c , and \x\<7q characterizes 
the amplitude of the trans-Planckian corrections to the fiducial spectrum. Martin & Ringeval 
(2004) report that the x 2 f° r suc h a model could give an improvement of 15 over the power 
law inflationary models for an additional 2 degrees of freedom with the first year WMAP 
data. With three years of data, many of the glitches and bites have disappeared, and the 
best fit trans-Planckian model of the form in equation (11) reduce the effective x 2 by only 4 
overall and by 5 in TT relative to power law inflation, a far less significant effect. 

The effect of the trans-Planckian corrections can be highly model dependent (See East- 
her et al. (2005a) and Easther et al. (2005b) for discussions). As an alternative, we consider 
forms that are more general as a way to look for oscillatory signals: 

A 2 ntTP (k) = A 2 ntPL (k) [1 + e TP cos 0(h)] (12) 

where 9 = v-j^ + or 9 = v\n + In these models, there are three new parameters: 
the amplitude, e^p, the frequency, v , and the phase, 0. 
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Assuming the ACDM model, we fit these three parameters to the data and find reduc- 
tions of 5 and 9.5 in the overall and TT xlff- As i n the Martin and Ringeval model, the 
improvements in the Xeff are driven by improvements in the fit around £ ~ 30 — 100 and 
the first peak. 



6. WMAP + Inflation 

The inflationary paradigm (Guth 1981; Sato 1981; Linde 1982; Albrecht & Steinhardt 
1982) explains the homogeneity, isotropy and flatness of the universe by positing an early 
epoch of accelerated expansion. This accelerated period of expansion also generated su- 
perhorizon fluctuations (Guth & Pi 1982; Starobinsky 1982; Mukhanov & Chibisov 1981; 
Hawking 1982; Bardeen et al. 1983). In the simplest inflationary models, these fluctua- 
tions are Gaussian distributed in amplitude with random phase and a nearly scale invariant 
spectrum of fluctuations. 

The detailed predictions of inflationary models depend on the properties of the inflaton 
potential (see Linde (2005) and Lyth & Riotto (1999) for recent reviews). Simple inflationary 
models predict that the slope of the primordial power spectrum, n s , differs from 1 and also 
predict the existence of a nearly scale-invariant spectrum of gravitational waves. In this 
section, we compare the simplest inflationary models to the WMAP three year data and 
to other cosmological data sets. Since we are constraining models with tensor modes, we 
also use the WMAP constraints on the amplitude of the B mode signal in the analysis. We 
characterize these models by seven basic parameters (the six basic parameters of the ACDM 
model plus one additional parameter, r, the ratio of the tensor to scalar power spectrum). 
Figure 14 shows the likelihood contours for the slope of the scalar fluctuations and the 
amplitude of the gravitational wave signal. 

The WMAP three year data place significant constraints on inflationary models. The 
strength of these constraints is apparent when we consider monomial models for the inflaton 
potential, V(<f>) oc <p a . These models (Lyth & Riotto 1999) predict 

r = loei ~ — 
N 

i-n. = 2e * + e ^^r (13) 

where N is the number of e-folds of inflation between the epoch when the horizon scale 
modes left the horizon and the end of inflation. Figure 14 compares the predictions of these 
monomial inflationary models to the data. For N = 60, A0 4 predicts r = 4/15, n s = 0.95, just 
at the outer edge of the 3a contour. For N = 50, A0 4 predicts r = 0.32, n s = 0.94. However, 
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Fig. 14. — Joint two-dimensional marginalized contours (68% and 95% confidence 
levels) for inflationary parameters (ro.oo2 5 n s ). We assume a power-law primor- 
dial power spectrum, dn s /d\nk = 0, as these models predict a negligible amount 
of running index, dn s /dhxk « — 10~ 3 . (Upper left) WMAP only. (Upper right) 
WMAP+SDSS. (Lower left) WMAP+2dFGRS. (Lower right) WMAP+CBI+VSA. 
The dashed and solid lines show the range of values predicted for monomial in- 
flaton models with 50 and 60 e-folds of inflation (equation (13)), respectively. 
The open and filled circles show the predictions of m 2 2 and A0 4 models for 50 
and 60 e-folds of inflation. The rectangle denotes the scale-invariant Harrison- 
Zel'dovich-Peebles (HZ) spectrum (n s = l,r = 0). Note that the current data 
prefers the m 2 (j) 2 model over both the HZ spectrum and the A</> 4 model by likeli- 
hood ratios greater than 12. (5x 2 > 5) 
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Table 7: Best Fit Inflationary Parameters (WMAP data only) 



Parameter 


ACDM + Tensor 


ACDM + Running ±Tensors 


Q b h 2 

h 

n s 

dn s /d\nk 

r 

T 
0-8 


0.0233 ±0.0010 

n 1 1Qc;+0.0094 
U.l±yo_o.0093 

0.787 ±0.052 
0.9841°;°! 
set to 
< 0.65 (95% CL) 
0.090 ±0.031 
0.702 ±0.062 


0.0219 ±0.0012 
0.128 ±0.011 
0.731 ±0.055 

1.16 ±0.10 
-0.085 ± 0.043 
< 1.1 (95% CL) 
0.108±g;gg 
0.712 ±0.056 



if we allow for non-minimal gravitational couplings, then the gravity wave predictions of 
these models are significantly reduced (Hwang & Noh 1998; Komatsu & Futamase 1999) 
and the models are consistent with the data. Alternatively, the m 2 4> 2 model is a good fit 
to the observations and its predicted level of gravitational waves, r ~ 0.13 — 0.16, is within 
range of upcoming experiments. 

If we restrict our analysis to a specific model and value of N, then the model specifies 
both r and n s . In Table 8 we show the best fit parameters for the m 2 <f) 2 potential with 
N = 60. Since n s and r are set to 0.9667 and 0.1333 for this model, there are only five 
free parameters (Q m h 2 , Q b h 2 , h, r, and A). The likelihood for the best fit m 2 <j) 2 is 12 times 
higher than the best fit "Harrison Zel'dovich-Peebles" model with n s = 1 and r = 0. To 
be conservative, this comparison uses the exact (res-4) likelihood code up to Z = 30 and the 
slightly lower point source correction discussed in Appendix A. If we cut the exact likelihood 
off at Z = 12, and use the original, higher point source value, the likelihood ratio is 15. 

Table 8: Best Fit Parameters for m 2 <p 2 with iV = 60 (WMAP data only) 



Parameter 


Mean 


Best Fit 


Q b h 2 


0.02218 ±0.00045 


0.0222 


Vt m h 


0.1232 ±0.0076 


0.125 


h 


0.741 ±0.029 


0.734 


T 


0.081 ±0.027 


0.0893 


0-8 


0.734 ± 0.046 


0.749 



In Peiris et al. (2003), we used the inflationary flow equations (Hoffman & Turner 2001; 
Kinney 2002) to explore the generic predictions of inflationary models. Here, we use the 
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slow-roll approximation to explore the implications of the data for inflationary models. The 
results of the third year analysis are consistent with the conclusions from the first year data: 
while the data rule out large regions of parameter space, there are also wide range of possible 
inflationary models consistent with our current data. One of the most intriguing features 
of Figure 14 is that the data now disfavors the exact Harrison-Zel'dovich-Peebles spectrum 
(n s = 1, r = 0). For power law inflationary models, this suggests a detectable level of gravity 
waves. There are, however, many inflationary models that predict a much smaller gravity 
wave amplitude. Alternative models, such as the ekpyrotic scenario (Khoury et al. 2001, 
2002) also predict an undetectable level of gravity waves. 

There are several different ways of expressing the constraints that the CMB data impose 
on inflationary models. These parameters can be directly related to observable quantities: 
n s — 1 = — 2ei — e 2 and r = 16ei. For the power law models, the WMAP bound on r implies 
that ei < 0.03 (95% C.L.). An alternative slow roll representation (see Liddle & Lyth (1992, 
1993)) uses 

e« 

Vv 

These parameters can be related directly to observables: r = 16e v and n s — 1 = — 6e v + 2r) v . 
Peiris et al. (2003) discusses various classes of models in slow roll parameter space. 

Models with a significant gravitational wave contributions, r ~ 0.3, make specific predic- 
tions for CMB and large-scale structure observations: (a) a modified temperature spectrum 
with more power at low multipoles; and (b) a lower amplitude of density fluctuations (for 
fixed CMB temperature fluctuations). For power law models, the strongest CMB constraints 
come from the shape of the temperature spectrum and the amplitude of density fluctuations. 
In order to fit the CMB data, models with higher r values require larger values of n s and 
land a lower amplitude of scalar fluctuations, which impacts large-scale structure predic- 
tions. Thus the strongest overall constraints on the tensor mode contribution comes from 
the combination of CMB and large-scale structure measurements (see Table 6). These strong 
limits rely on our assumption of a power law spectral index. If we allow for a running index, 
then models with a large tensor components are consistent with the data. 



M 2 Pl 



XL 
v 

Yl 

V 



(14) 
(15) 
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Table 9: Constraints on r, Ratio of Amplitude of Tensor Fluctuations to Scalar Fluctuations 
(at k = 0.002 Mpc' 1 ) 



Data Set 


r (no running) 


r (with running) 


WMAP 
WMAP+BOOM+ACBAR 
WMAP+CBI+VSA 
WMAP+2df 
WMAP+SDSS 


< 0.65 (95% CL) 
<0.68 (95% CL) 
<0.62 (95% CL) 
<0.38 (95% CL) 
<0.30 (95% CL) 


<1.1 (95% CL) 
<1.1 (95% CL) 
<1.1 (95% CL) 
<0.64 (95% CL) 
<0.38 (95% CL) 



7. Constraining the Composition of the Universe 

7.1. Dark Energy Properties 

Over the past two decades, there has been growing evidence for the existence of dark 
energy (Peebles 1984; Turner et al. 1984; Ostriker & Steinhardt 1995; Dunlop et al. 1996; 
Krauss & Turner 1999; Bahcall et al. 1999). By measuring both the acceleration (Riess et al. 
1998; Perlmutter et al. 1999) and deceleration (Riess et al. 2004) of the universe, supernova 
observations provide the most direct evidence for the existence of dark energy. 

The nature of this dark energy is a mystery. From a field theoretic perspective the 
most natural explanation for this would be the presence of a residual vacuum energy density 
or cosmological constant, A, (Carroll et al. 1992; Peebles & Ratra 2003). However, there 
are well-known fine-tuning and coincidence problems in trying to explain the 120 orders- 
of-magnitude discrepancy between the expected "natural" Planck-scale energy density of 
a cosmological constant and the observed dark energy density. These problems motivate 
a wide range of alternative explanations for the observations including the presence of an 
extra matter candidate: for example a dynamical, scalar "quintessence" field (Peebles & 
Ratra 1988; Wetterich 1988; Zlatev et al. 1999), minimally coupled (Caldwell et al. 1998; 
Ferreira & Joyce 1998) or non-minimally coupled to gravity (Amendola 1999) or other matter 
(Bean & Magueijo 2001). In this latter case, the measured acceleration is due to interactions 
in the matter bulk. Another alternative is that modifications to gravity (e.g., Deflayet (2001); 
Deffayet et al. (2002)) are responsible for the observed anomalies. 

The dark energy has two distinct cosmological effects: (1) through the Friedman equa- 
tion, it alters the evolution of H (z) and (2) through the perturbation equations, it alters the 
evolution of D(z), the growth rate of structure. The supernova data measure only the lumi- 
nosity distance, which depends on H(z), while the large scale structure data are sensitive to 
both H(z) and D(z). 
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Fig. 15. — Constraints on w, the equation of state of dark energy, in a flat uni- 
verse model based on the combination of WMAP data and other astronomical 
data. We assume that w is independent of time, and ignore density or pressure 
fluctuations in dark energy. In all of the figures, WMAP only constraints are 
shown in blue and WMAP + astronomical data set in red. The contours show 
the joint 2-d marginalized contours (68% and 95% confidence levels) for Q m and 
w. ( Upper left) WMAP only and WMAP + SDSS. ( Upper right) WMAP only and 
WMAP + 2dFGRS. (Lower left) WMAP only and WMAP+SN(HST/GOODS). 
(Lower right) WMAP only and WMAP+SN(SNLS). In the absence of dark en- 
ergy fluctuations, the excessive amount of ISW effect at i < 10 places significant 
constraints on models with w < —1. 



While the presence of dark energy impacts the CMB primarily through the distance to 
the surface of last scatter, the dark energy clustering properties also affect the CMB prop- 
erties. The dark energy response to gravitational perturbations depends upon its isotropic 
and anisotropic sound-speeds (Hu 1998; Bucher & Spergel 1999). This affects the CMB fluc- 
tuations through the ISW effect. If the dark energy can cluster, then it produces a smaller 
ISW effect and does not enhance the power spectrum at large angular scales. These effects 
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Table 10: Constraints on w in Flat Cosmologies With Different Assumption About Dark 
Energy Clustering 



Data Set 


with perturbations 


no perturbations 


WMAP + SDSS 
WMAP + 2dFGRS 
WMAP + SNGold 
WMAP + SNLS 
CMB+ LSS+ SN 


-0.75^ 6 7 

-0.8918:12 

-0 qiq+ - 081 

u.y±»_ QgQ 

-0 q67 +0 ' 073 
u.yu i _ .072 

-o q2fi +0 - 054 

U.»z,u_ 053 


-0.69 ±0.16 

n O77+0.097 
-U.8C '_ .098 

-0 q42 +0 ' 075 
u.y^_ 073 

-0 q84+0- 070 
-0.915 ±0.051 



are most dramatic for models with w < — 1, as dark energy effects in these models turn on 
suddenly at late times and significantly enhance the quadrupole. This can be understood in 
terms of the constraints imposed by the shape of the angular power spectrum: if we assume 
that the dark energy properties can be described by a constant value of w, then fixed peak 
position and heights (which determine VL m h 2 ) confine our models to a narrow valley in the 
(Q m ,w) likelihood surface as shown in Figure 15 and 16. The figures show that the 3 year 
data enable a more accurate determination of Q m h 2 which narrows the width of the degen- 
eracy valley. The pair of figures show that CMB data can place strong limits on models with 
w < — 1 and non-clustering dark energy. On the other hand, if the dark energy is a matter 
component that can cluster, even meagerly, as is the case for scalar field theories where c^—l, 
then this clustering counters the suppression of perturbation growth during the accelerative 
epoch and the quadrupole's magnitude is reduced. This lessens the discriminating power of 
the quadrupole for measuring w: while CMB data rule out the w << — 1 region in Figure 
15, it does not constrain models in the same region in Figure 16. 

It's interesting to note that if we relax the assumption of spatial flatness, the data still 
prefers models with w close to —1 (see Figure 17). In this analysis, we assume that the dark 
energy clusters. 
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Fig. 16. — Constraints on w, the equation of state of dark energy, in a fiat universe, 
Model M6 in Table 3, based on the combination of WMAP data and other 
astronomical data. We assume that w is independent of time, but include density 
and pressure fluctuations in dark energy with the speed of sound in the comoving 
frame equal to the speed of light, c 2 s = \. In all of the figures, WMAP data only 
constraints are shown in black solid lines and WMAP + astronomical data set 
in red. The contours show the joint 2-d marginalized contours (68% and 95% 
confidence levels) for Vt m and w. ( Upper left) WMAP only and WMAP + SDSS. 
( Upper right) WMAP only and WMAP + 2dFGRS. (Lower left) WMAP only and 
WMAP+SNgold. (Lower right) WMAP only and WMAP+SNLS. In the presence 
of dark energy fluctuations, the ISW effect at i < 10 is nearly canceled by dark 
energy fluctuations and thus the WMAP data alone do not place significant 
constraints on w. The WMAP only models are shown assuming no prior on Ho 
(dashed black) and with an assumed prior of H < 100km/ s/Mpc (full black lines). 
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Fig. 17. — Constraints on a non-flat universe with quintessence-like dark energy 
with constant w (Model M10 in Table 3). The contours show the 2-d marginalized 
contours for w and £lk based on the the CMB+2dFGRS+SDSS+supernova data 
sets. This figure shows that with the full combination of data sets, there are 
already strong limits on w without the need to assume a flat universe prior. 
The marginalized best fit values for the equation of state and curvature are 
w = -1.08 ± 0.12 and Vl k = -0.026l°:oi5 at the 68% confidence level. 
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Fig. 18. — Constraints on a flat universe with quintessence-like dark energy and 
non-relativistic neutrinos. The contours show the 2-d marginalized constraints 
on r the mass of non-relativistic neutrinos, m u , and the dark energy equation of 
state, w, assumed constant, based on the the CMB+2dFGRS+SDSS+supernova 
data sets. In this analysis, we assume that the dark energy clusters. With 
the combination of CMB+2dFGRS+SDSS+supernova data sets, there is not a 
strong degeneracy between neutrino and dark energy properties. Even in this 
more general model, we still have an interesting constraint on the neutrino mass 
and equation of state: < 1.0 eV (95% CL) and w = -1.07 ±0.12 (68% CL). 

This suggests that the astronomical dark energy and neutrino limits are robust 
even when we start to consider more baroque models. 
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7.2. Neutrino Properties 

7.2.1. Neutrino Mass 

Both atmospheric neutrino and solar neutrino experiments show that neutrinos are mas- 
sive and that there is significant mixing between the various neutrino interaction eigenstates 
(see Mohapatra et al. (2005) for a recent review). These experiments measure the difference 
between the square of the neutrino masses, m 2 u . — m% , rather than the mass of individual 
neutrino mass eigenstates. Cosmological measurements nicely complement these measure- 
ments by constraining £\ m Vi . Since light massive neutrinos do not cluster as effectively as 
cold dark matter, the neutrino mass has a direct impact on the amplitude and shape of the 
matter power spectrum (Bond et al. 1980; Bond & Szalay 1983; Ma 1996; Hu et al. 1998). 
The presence of a significant neutrino component lowers the amplitude of matter fluctuations 
on small scales, a by roughly a factor proportional to m„), where m v is the total mass 
summed over neutrino species, rather than the mass of individual neutrino species. The 
current large-scale structure data restrict A In erg < 0.2, but they are not sensitive enough to 
resolve the free-streaming scale of individual neutrino species (Takada et al. 2005). 

Using a combination of the first year WMAP data, small-scale CMB and large-scale 
structure data, Spergel et al. (2003) placed an upper limit on X]« m ^ < 0-7 e V- While 
this limit does not depend on the Lyman a data, it is sensitive to galaxy bias measurements 
(which normalizes the large-scale structure data) and to the addition of small scale CMB data 
(which improves the measurement of cosmological parameters). Over the past year, several 
groups obtained comparable (but slightly different) limits (Hannestad 2003; Pierpaoli 2003; 
Elgar0y & Lahav 2003). The differences are due to including (or removing) external data 
sets and priors or to adding additional cosmological parameters. 

The limits on neutrino masses from WMAP data alone is now very close to limits based 
on combined CMB data sets. Ichikawa et al. (2005) used the CMB data alone to place a limit 
on the neutrino mass of ^2 m v < 2.0 eV. Using WMAP data alone, we now find Yl m u < 1-8 
eV. 

Since the presence of massive neutrinos slows the growth of small scale structure, the 
combination of CMB and large-scale structure data constrain the neutrino mass. Figure 
19 shows the likelihood function as a function of neutrino mass and amplitude of mass 
fluctuations in the local universe, a%. The 95% confidence limits on neutrino mass are 
given in Table 11. The combination of WMAP with SDSS and WMAP with 2dFGRS data 
constrain a s at roughly the same level, 20% at the 95% confidence level. This constraint 
yields comparable limits on the neutrino mass. While the WMAP data have improved, the 
error bars on <jg have not significantly changed from the limits obtained from WMAPext + 
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2dFGRS, thus, the limit on neutrino mass is quite close to the first year limit. Note that in 
the first year analysis, we used the (Verde et al. 2002) measurement of bias for the 2dFGRS 
preliminary data as there had not been an equivalent analysis done for the full 2dFGRS 
data set. As discussed in §4.1.4, we now marginalize over the 2dFGRS bias and use the bias 
measurements of (Seljak et al. 2005b) for SDSS. 

If the constraints on amplitude are reliable, then small scale matter power spectrum 
structure data can significantly improve these neutrino constraints. Goobar et al. (2006) 
have recently completed a CMB + Lyman a study and place a limit of s ^ J m v < 0.30eV 
(95% C.L.). Similarly, cluster-based measurements of as and lensing-based measurements of 
as have the potential to tighten the constraint on m v . 

0.9 r — ' — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — ; 
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Fig. 19. — Joint two-dimensional marginalized contours (68% and 95% confidence 
levels) on (a 8 ,m u ) for WMAP only, Model M7 in Table 3, and WMAP+SDSS. By 
measuring the growth rate of structure from z = 1088 to z ~ 0, these observations 
constrain the contribution of non-relativistic neutrinos to the energy density of 
the universe. 
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Table 11: Constraints on Neutrino Properties 



Data Set 


Y,m v (95% limit for N v = 3.04) 


N v 


WMAP 
WMAP + SDSS 
WMAP + 2dFGRS 
CMB + LSS +SN 


1.8 eV (95% CL) 
1.3 eV (95% CL) 
0.88 eV (95% CL) 
0.66 eV (95% CL) 


7 1+ 4 - 1 

'• 1 -3.5 

2.7 ± 1.4 
3.3 ± 1.7 



7.2.2. Number of Relativistic Species 

If there are other light stable neutral particles (besides the three light neutrinos and 
the photon), then these particles will affect the CMB angular power spectrum and the 
evolution of large-scale structure. Because of the details of freeze-out at electron-positron 
annihilation (Gnedin & Gnedin 1998; Dolgov, Hansen & Semikov 1999), QED corrections 
at finite temperature (Mangano et al. 2002) and non-trivial terms in the neutrino mass 
matrix (Mangano et al. 2005), the effective number of neutrino species is slightly greater 
than 3. Any light particle that does not couple to electrons, ions and photons will act as 
an additional relativistic species. For neutrinos, we can compute their effects accurately as 
their temperature is (4/1 1) 1 / 3 of the CMB temperature. For other relativistic species, the 
limit on N* ff - 3.04 can be converted into a limit on their abundance by scaling by the 
temperature. 

The shape of the CMB angular power spectrum is sensitive to the epoch of mat- 
ter/radiation equality. If we increase N v , the effective number of neutrino species, then 
we will need to also increase the cold dark matter density, Q c h 2 , and slowly change other 
parameters to remain consistent with the WMAP data (Bowen et al. 2002). In addition, the 
presence of these additional neutrino species alters the damping tail and leaves a distinctive 
signature on the high £ angular power spectrum (Bashinsky & Seljak 2004) and on the small 
scale matter power spectrum. 

The high matter density also alters the growth rate of structure, thus, the combination 
of large-scale structure and CMB data constrains the existence of any new light relativistic 
species. These limits constrain both the existence of new particles and the interaction prop- 
erties of the neutrino (Bowen et al. 2002; Hall & Oliver 2004). Hannestad (2001) used the 
pre- WMAP CMB and large-scale structure data to place an upper limit of N v < 17. After 
the release of the first year WMAP data, several authors (Hannestad 2003; Pierpaoli 2003; 
Barger et al. 2003; Crotty et al. 2003; Elgar0y & Lahav 2003; Barger et al. 2004; Hannestad 
2005) used the combination of WMAP, 2dFGRS and various external data to reduce this 
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limit by a factor of 2-3. Table 11 shows the maximum likelihood estimate of the number of 
neutrino species for different data set combinations using the new WMAP data. The SDSS 
and 2dFGRS data differ in the shapes of the two measured power spectra: this difference 
leads to the disagreement in their best fit values for N^ff. 



7.3. Non-Flat Universe 



The WMAP observations place significant constraints on the geometry of the universe 
through the positions of the acoustic peaks. The sound horizon size, r s , serves as a very useful 
ruler for measuring the distance to the surface of last scatter. For power law open universe 
models, 6 A = 0.535+°;°^ °. Figure 21 shows that this constraint confines the likelihood 
function to a narrow degeneracy surface in (fl m ,fl\). This degeneracy line is well fit by 




Fig. 20. — Range of non-flat cosmological models consistent with the WMAP data 
only. The models in the figure are all power-law CDM models with dark energy 
and dark matter, but without the constraint that VL m + Vt^ = 1 (model M10 in 
Table 3). The different colors correspond to values of the Hubble constant as 
indicated in the figure. While models with Q\ = are not disfavored by the 
WMAP data only (Axl ff = 0; Model M4 in Table 3), the combination of WMAP 
data plus measurements of the Hubble constant strongly constrain the geometry 
and composition of the universe within the framework of these models. The 
dashed line shows an approximation to the degeneracy track: Qk = —0.3040 + 
0.4067^a- Note that for these open universe models, we assume a flat prior on 
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Fig. 21. — Joint two-dimensional marginalized contours (68% and 95%) for mat- 
ter density, Q m , and vacuum energy density, for power-law CDM models 
with dark energy and dark matter, but without the constraint that Q m + ^a — 1 
(model M10 in Table 3). The panels show various combinations of WMAP and 
other data sets. While models with Vt m = 0.415 and = 0.630 are a better fit 
to the WMAP three year data alone than the flat model, the combination of 
WMAP three year data and other astronomical data favors nearly flat cosmolo- 
gies. ( Upper left) WMAP+HST key project measurement of H . ( Upper right) 
WMAP+SDSS LRG measurement of the angular diameter distance to z = 0.35. 
(Middle left) WMAP+SNLS data. (Middle right) WMAP+SNGold. (Lower left) 
WMAP+2dFGRS. (Lower right) WMAP+SDSS. 
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Q K = —0.3040 + 0.4067f2A. However, the CMB data alone does not distinguish between 
models along the valley: it is consistent with both flat models and models with = 0. If 
we allow for a large SZ signal, then the WMAP data alone favors a model with — —0.04; 
however, this model is not consistent with other astronomical data. 

The combination of WMAP data and other astronomical data places strong constraints 
on the geometry of the universe (see Table 12): 

• The angular scale of the baryon acoustic oscillation (BAO) peak in the SDSS LRG 
sample (Eisenstein et al. 2005) measures the distance to z — 0.35. The combination of 
the BAO and CMB observations strongly constrain the geometry of the universe. The 
position of the peak in the galaxy spectrum in the SDSS and 2dFGRS surveys provide 
local measurements of the angular diameter distance. 

• Figure 20 shows that the Hubble constant varies along this line, so that the HST key 
project constraint on the Hubble constant leads to a strong bound on the curvature. 

• SNe observations measure the luminosity distance to z ~ 1. The combination of SNe 
data and CMB data also favors a nearly flat universe. 

The strong limits quoted in Table 12 rely on our assumption that the dark energy has 
the equation of state, w — — 1. In section 7.1, we discussed relaxing this assumption and 
assuming that w is a constant. Figure 15 shows that by using the combination of CMB, 
large-scale structure and supernova data, we can simultaneously constrain both Q k and w. 
This figure confirms that our minimal model, Q k = and w = — 1 is consistent with the 
current data. 

Table 12: Joint Data Set Constraints on Geometry and Vacuum Energy 



Data Set 






WMAP + h = 0.72 ± 0.08 
WMAP + SDSS 
WMAP + 2dFGRS 
WMAP + SDSS LRG 
WMAP + SNLS 
WMAP + SNGold 


-0.014 ±0.017 

-o no c i3 +0 0068 

U.UUOO_ 0060 

-0 nnQS+ - 0098 
u.uuyo_o.oo92 

-0.012 ±0.010 
-0.011 ±0.012 
-0.023 ±0.014 


0.716 ±0.055 

0.707 ±0.041 
n 74C+0.025 

0.728 ±0.021 
0.738 ±0.030 
0.700 ±0.031 
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8. Are CMB Fluctuations Gaussian? 

The detection of primordial non-Gaussian fluctuations in the CMB would have a pro- 
found impact on our understanding of the physics of the early universe. While the simplest 
inflationary models predict only mild non-Gaussianities that should be undetectable in the 
WMAP data, there are a wide range of plausible mechanisms for generating significant and 
detectable non-Gaussian fluctuations (Bartolo et al. (2004a) for a recent review). There are 
a number of plausible extensions of the standard inflationary model (Lyth et al. 2003; Dvali 
et al. 2004; Bartolo et al. 2004b) or alternative early universe models (Arkani-Hamed et al. 
2004; Alishahiha et al. 2004) that predict skewed primordial fluctuations at a level detectable 
by WMAP. 

There are other cosmological mechanisms for generating non-Gaussianity. The smallness 
of the CMB quadrupole seen by both WMAP and COBE has stimulated interest in the 
possibility that the universe may be finite (Luminet et al. 2003; Aurich et al. 2005). If 
the universe were finite and had a size comparable to horizon size today, then the CMB 
fluctuations would be non-Gaussian (Cornish et al. 1996; Levin et al. 1997; Bond et al. 2000; 
Inoue et al. 2000). While analysis of the first year data did not find any evidence for a finite 
universe (Phillips & Kogut 2006; Cornish et al. 2004), these searches were non-exhaustive so 
the data rule out most but not all small universes. 

Using an analysis of Minkowski functionals, Komatsu et al. (2003) did not find evidence 
for statistically isotropic but non-Gaussian fluctuations in the first year sky maps . The 
Colley & Gott (2003) reanalysis of the maps confirmed the conclusion that there was no 
evidence of non-Gaussianity. Eriksen et al. (2004a) measured the Minkowski functionals 
and the Length of the Skeleton for the first year maps on 11 different smoothing scales. 
While they found no evidence for deviations from non-Gaussianity using the Minkowski 
area, Minkowski length and the length of the skelton, they did find a intriguingly high \ 2 
for the genus statistic. 

For a broad class of theories, we can parameterize the effects of non-linear physics by 
a simple coupling term that couples a Gaussian random field, ip, to the Bardeen curvature 
potential, $: 

$(£) = + W 2 (x) (16) 

Simple inflationary models based on a single slowly-rolling scalar field with the canoni- 
cal kinetic Lagrangian predict |/jvx| < 1 (Maldacena 2003; Bartolo et al. 2004a); how- 
ever, curvaton inflation (Lyth et al. 2003), ghost inflation (Arkani-Hamed et al. 2004), and 
Dirac-Born-Infeld (DBI) inflation models (Alishahiha et al. 2004) can generate much larger 
non-Gaussianity, |/jvl| ~ 100. Using the WMAP first year data, Komatsu et al. (2003) con- 
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strained —54 < f NL < 134 at the 95% confidence level. Several different groups (Gaztanaga 
& Wagg 2003; Mukherjee & Wang 2003; Cabella et al. 2004; Phillips & Kogut 2006; Crem- 
inelli et al. 2005) have applied alternative techniques to measure f^L from the maps and 
have similar limits on /jvl- Babich et al. (2004) note that these limits are sensitive to the 
physics that generated the non-Gaussianity as different mechanisms predict different forms 
for the bispectrum. 

Since the release of the WMAP data, several groups have claimed detections of signif- 
icant non-Gaussianities (Tegmark et al. 2003; Eriksen et al. 2004b; Copi et al. 2003; Vielva 
et al. 2004; Hansen et al. 2004; Park 2004; Larson & Wandelt 2004; Cruz et al. 2005). 
Almost all of these claims imply that the CMB fluctuations are not stationary and claim 
a preferred direction or orientation in the data. Hajian et al. (2005) argue that these de- 
tections are based on a posteriori selection of preferred directions and do not find evidence 
for preferred axes or directions. Because of the potential revolutionary significance of these 
detections, they must be treated with some caution. Galactic foregrounds are non-Gaussian 
and anisotropic, and even low level contamination in the maps can produce detectable non- 
Gaussianities (Chiang et al. 2003; Naselsky et al. 2005), but have only minimal effects on 
the angular power spectrum (Hinshaw et al. 2003). Similarly, point sources can be a source 
of non-Gaussianity at small angular scales (Eriksen et al. 2004b). Because of the WMAP 
scan pattern, the variance in the noise in the maps is spatially variable. There is significant 
1/ / noise in several of the Difference Assemblies (DAs) (particularly W4) — which leads to 
anisotropies in the two-point function of the noise. Finally, most of the claimed detections 
of significant non-Gaussianities are based on a posteriori statistics. Many of the claimed 
detections of non-Gaussianity can be tested with the three year WMAP data (available at 
lambda.gsfc.nasa.gov). Future tests should use the simulated noise maps, Monte Carlo sim- 
ulations and the difference maps (year 1 — year 2, year 2 — year 3, etc.) to confirm that the 
tests are not sensitive to noise statistics and the multi-frequency data to confirm that any 
claimed non-Gaussianity has a thermal spectrum. Claims of non-Gaussianity incorporating 
data close to the galactic plane (within the Kp2 cut) should be treated with caution, as the 
foreground corrections near the plane are large and uncertain. 

The following subsections describe a number of statistical tests designed to search for 
non-Gaussianities in the microwave sky. All of these analyses use three year maps cleaned 
with the KKaHaDust templates (Hinshaw et al. 2006). We refer to these maps as the 
"template-cleaned maps" . In the first subsection, we show that the probability distribution 
function of the cleaned CMB maps is consistent with Gaussianity. In the second subsec- 
tion, we show that the Minkowski functionals are consistent with expectations for Gaussian 
fluctuations. Next, we compute the bispectrum of the cleaned maps. The final subsection 
describes measurements of the four point function. 
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8.1. One Point Distribution Function 

One of the simplest tests of non-Gaussianity is a measurement of the one point prob- 
ability function. However, because the detector noise in the map is inhomogeneous (higher 
in the ecliptic plane and lower near the poles), this test is non-trivial. We account for the 
spatial variations in noise by computing a variance-normalized temperature for each pixel in 
a given map: 

Ui = Ti (17) 

V °~noise/ N obs + &CMB 

where Tj is the measured temperature signal, the detector noise depends on the number of 
observations of a given pixel, N j, s . Here, we apply the analysis to template-cleaned maps 
outside the Kp2 skycut. For this analysis, we compute <J no i se , the noise per observation, from 
the year 1 — year 2 difference maps and fit cr CMB , the CMB signal, to the sum of the year 
one and year two maps. With N side = 1024, the computed a noise value is within 0.5% of the 
value of (T estimated from the time series (Jarosik et al. 2006). As we lower the resolution, 
the value of a no ise slowly increases with increasing pixel size. For W4, the channel with the 
large 1/ / noise, this change is most dramatic; the value of a at resolution N side = 32 is 6% 
higher than the value computed for N side = 1024. 

Figures 22 and 23 shows the one-point distribution function of the cleaned sky maps as 
a function of resolution. Most of the distributions appear to be rather well fit by a Gaussian. 
This result is consistent with that from the area of hot and cold spots (one of the Minkowski 
functionals), which measures the cumulative one point probability function. However, with 
N S ide = 16, there is a slight excess at negative T in the first panel of Figure 23. This is due 
to the cold spot near I = 209 and b = —57 discussed by (Vielva et al. 2004) in the context of 
the first year data. Cruz et al. (2006) argue that there is a 1.85% probability of seeing the 
level of skewness and kurtosis associated with this cold spot. While this probability would 
increase if they allowed for variable smoothing scales, the results are intriguing and merit 
further investigation. 
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Fig. 22. — Normalized one point distribution function of temperature anisotropy, 
defined in equation (17), for the template-cleaned Q (left), V (middle) and W 
(right) band maps outside the Kp2 cut. The sky maps have been degraded to 
N S ide — 256 for this figure. The red line shows the Gaussian distribution, which 
is an excellent fit to the one point distribution function. 




Fig. 23. — Normalized one point distribution function of temperature anisotropy, 
defined in equation (17), for the template-corrected V band data maps outside 
the KpO cut. The sky maps have been degraded to N side = 16(left), 64(middle) 
and 256(right) for this figure. The red line shows the best fit Gaussian, which is 
an excellent fit to the one point distribution function. 
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8.2. Size and Shape of Hot and Cold Spots 

Minkowski functionals (Minkowski 1903; Gott et al. 1990; Mecke, Buchert & Wagner 
1994; Schmalzing & Buchert 1997; Schmalzing & Gorski 1998; Winitzki & Kosowsky 1998) 
measure the contour length, area, and number of hot and cold spots. Following the approach 
used in the first year analysis, we compute the Minkowski functionals as a function of tem- 
perature threshold, v = AT /a, where a is the standard deviation of the map. For a two 
dimensional map, we measure three Minkowski functionals, the genus, G(u), of the maps, 
the contour length, C(u) and the area within the contours, A(v). 

We compare the measured values of the Minkowski functionals to their expected am- 
plitude for a Gaussian sky. We simulate a series of maps based on the best fit parameters 
for ACDM and the WMAP noise patterns. For the analysis, we use the template-cleaned 
V+W maps outside the KpO sky mask region. Following the approach used in Komatsu 
et al. (2003), we compute the Minkowski functionals at 15 thresholds from — 3.5a to +3.5cx 
and compare each functional to the simulations using a goodness of fit statistic, 

X 2 = ^2 [ F WMAP - ( F sim)] Ul ^uiu 2 [ F WMAP ~ { F lim)] V2 ( 18 ) 

where Fwmap * s the Minkowski functional computed from the WMAP data, F % sim is the 
Minkowski functional computed from the simulated data, and £„ 1I/2 is the bin-to-bin co- 
variance from the simulations. Figure 24 shows the Minkowski functionals as a function of 
threshold for a map with N side = 128 (28 ' pixels). These pixels are small enough to resolve 
the acoustic spots, but not so small as to be noise dominated. The figure shows that the con- 
tour length, area, and number of spots is consistent with expectations for a Gaussian theory. 
Table 13 lists the probability of measuring the observed values of the Minkowski function- 
als as a function of pixel size. At all resolutions, the maps are consistent with Gaussian 
simulations. 

We have also simulated non-Gaussian sky with non-Gaussian signals generated according 
to equation (16). By comparing these simulations to the data, we can constrain f^L = 7±66 
at the 68% confidence level, consistent with the bispectrum measurement (§8.3). 

8.3. Bispectrum 

The bispectrum is sensitive to both primordial non-Gaussianity and various sources 
of secondary anisotropy (Spergel & Goldberg 1999; Goldberg & Spergel 1999; Komatsu & 
Spergel 2001). Here, we use the WMAP 3 year data to constrain the amplitude of primordial 



-64- 






Fig. 24. — The WMAP data are in excellent agreement with the Gaussian sim- 
ulations based on the analysis of the Minkowski functionals for the three year 
WMAP data outside the KpO cut. The filled circles in the left panel shows the 
values for the data at A^e = 128 (28' pixels). The gray band shows the 68% 
confidence interval for the Gaussian Monte Carlo simulations. The right panels 
show the residuals between the mean of the Gaussian simulations and the WMAP 
data. Note that the residuals are highly correlated from bin to bin, so the x 2 are 
consistent with the Gaussian simulations. 
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Table 13: \ 2 for Minkowski Functionals for 15 thresholds for the template-cleaned VW 
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Minkowski 
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20.9 
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19.2 


15 
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14.0 


15 


15.3 


0.54 
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45 
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0.31 
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15 
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15 
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0.93 
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17.3 
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15 
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32 
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15 
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32 
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43.8 


45 


49.1 


0.61 
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28.2 


15 


15.8 
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16 
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15 
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0.29 


16 
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14.1 


15 


15.6 
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16 
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84.6 


45 


49.4 


0.03 


8 
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10.8 


15 


15.5 


0.62 


8 


Contour 


24.3 


15 


16.0 


0.09 


8 


Spot Area 


28.8 


15 


15.0 


0.05 


8 


Combined 


100.5 


45 


49.0 


0.03 



non-Gaussianity and to detect the amplitude of the point source signal in the cleaned Q, V 
and W band maps. 

The amplitude of the primordial non-Gaussian signal can be found by computing a cubic 
statistic on the map (Komatsu et al. 2005): 



5 'primordial = j^- J ^r 2 dr J ^ A(r , n) B 2 (r , n) (19) 



where f s k y is the fraction of the sky used in the analysis, B(r, n) is a Weiner filtered map 
of the primordial fluctuations and A is optimized to detect the form of the non-linearities. 
The amplitude of S pr i mor di a i can be related directly to /atl- Here, we use A and B as defined 
in Komatsu et al. (2005). While we used i max = 265 for the first-year analysis, we use 
^max = 350 for the present analysis, as noise is significantly lower with three years of data. 
The error on f NL begins to increase at t\ max > 350 due to the presence of inhomogeneous 
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noise. Note that Creminelli et al. (2005) argue that the optimal estimator for S pr i mor di a i 
should include a term that is linear in temperature anisotropy as well as a cubic term that 
we already have in equation (19). They claim that their estimator could reduce the error on 
I'nl by about 20%. While their result is attractive, we shall not include the linear term in 
our analysis, as their estimator has not been tested against non-Gaussian simulations and 
thus it is not yet clear if it is unbiased. 

Point sources are an expected cause of non-Gaussianity. Because point sources are 
not very correlated on the angular scales probed by WMAP, the point sources make a 
constant contribution to the bispectrum, b src . Komatsu et al. (2005) develops a cubic statistic 
approach for computing b src : 



where m 3 = (An) 1 J d 2 hM 3 (h), M(n) = [(Tq MB + N(n)}~ 1 , and D(n) is a match filter 
optimized for point source detection: 



where be is a beam transfer function and Ce = C^ mb b 2 + iVj. We weight the temperature 
maps by M(n) before we calculate a\ m . We use i m ax = 1024 for calculating D(n). (See § 3.2 
of Komatsu et al. (2003) for details of the weighting method.) Given the uncertainties in 
the source cut-off and the luminosity function, the values for b src in Table 14 are consistent 
with the values of A ps in Hinshaw et al. (2006). 

Table 14 lists the measured amplitude of the non-Gaussian signals in the 3 year maps. 
The values are computed for template-cleaned Q, V and W band maps. With three years 
of data, the limits on primordial non-Gaussianity have improved from —58 < /tvl < 137 
to —54 < f NL < 114 at the 95% confidence level. The improvement in limit on f NL is 
roughly consistent with the expectation that in the signal-dominated regime, Af^L oc l^ ax 
(Komatsu & Spergel 2001). The level of point source non-Gaussianity in the 3 year maps 
is lower than in the first year maps. This drop is due to the more sensitive 3 year masks 
removing additional sources. 
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8.4. Trispectrum 



Motivated by claims that there are large scale variations in the amplitude of fluctua- 
tions, we consider a non-Gaussian model that generates a non-trivial four point function for 
the curvature (and temperature) fluctuations, but does not produce a three-point function. 
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Fig. 25. — Constraints on the amplitude of the four point function. The measured 
amplitude of the four point function (expressed in terms of a non-Gaussian 
amplitude defined in equation (23)) is compared to the same statistic computed 
for simulated Gaussian random fields. The yellow line shows the results for Q, 
V and W bands and the red histogram shows the distribution of the results of 
the Monte-Carlo realizations. Note that in both the simulations and the data A 
is greater than due to the inhomogeneous noise. The excess in Q is may be 
due to point source contamination. 



This model describes a cosmology where the value of one field modulates the amplitude of 
fluctuations in a second field: 

$(£) = 0(f) [1 + W@)] ( 22 ) 

where <fi and ip are Gaussian random fields and $ is the Bardeen curvature potential. The 
presence of such a term would generate variations in the amplitude of fluctuations across the 
sky. 

Appendix B derives an estimator for the amplitude of the non-Gaussian term, g 2 NL \ip\ 2 - 
This estimator is based on approximating the CMB fluctuations as arising from an infinitely 



Table 14: Amplitude of Non-Gaussianity 





/nl 


^src 

[10" 5 /iK 3 sr 2 ] 
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41 ±55 


4.8 ±2.0 


V 


25 ±50 


0.12 ±0.52 


w 


11 ±50 


-0.21 ±0.34 


v+w 


18 ±46 


0.25 ±0.26 


Q+V+W 


30 ±42 


0.73 ±0.36 
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thin surface of last scatter. We measure the amplitude of the four point function by com- 
puting 

Q = ECtfV 2 ?/ - ( 23 ) 

i 

where is a smoothed map (e.g., an N sil i e = 128 map), Tj is an unsmoothed map, and iVj 
is the expected value of Ty V 2 Tj for a map without any signal. 

Figure 25 shows measurements of Q from the Q, V and W band data. The V and W 
band data show any evidence for a non-trivial four point function, while Q band data may 
show the contamination from point sources. At the S/N level of the 3 year data, there are no 
significant cosmological or systematic effects modulating the amplitude of the fluctuations 
as a function of scale. 



9. Conclusions 

The standard model of cosmology has survived another rigorous set of tests. The errors 
on the WMAP data at large I are now three times smaller and there has been significant 
improvements in other cosmological measurements. Despite the overwhelming force of the 
data, the model continues to thrive. This was the basic result of (Spergel et al. 2003) and 
was reinforced by subsequent analyses of the first year WMAP data with the SDSS (Tegmark 
et al. 2004a), and analysis of first year WMAP plus the final 2dFGRS survey (Sanchez et 
al. 2006). After the analyses in this paper were completed, a larger SDSS LRG sample was 
released. An analysis of the three year WMAP data combined with the new larger sample 
reached the same basic conclusions as this paper: the ACDM model is consistent with both 
CMB and large-scale structure measurements and its basic parameters are now reasonably 
well constrained (Tegmark et al. 2006). 

The data are so constraining that there is little room for significant modifications of 
the basic ACDM model. The combination of WMAP measurements and other astronomical 
measurements place significant limits on the geometry of the universe, the nature of dark 
energy, and even neutrino properties. While allowing for a running spectral index slightly 
improves the fit to the WMAP data, the improvement in the fit is not significant enough to 
require a new parameter. 

Cosmology requires new physics beyond the standard model of particle physics: dark 
matter, dark energy and a mechanism to generate primordial fluctuations. The WMAP data 
provides insights into all three of these fundamental problems: 

• The clear detection of the predicted acoustic peak structure implies that the dark 
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matter is non-baryonic. 

• The WMAP data are consistent with a nearly flat universe in which the dark energy has 
an equation of state close to that of a cosmological constant, w — —1. The combina- 
tion of WMAP data with measurements of the Hubble Constant, baryon oscillations, 
supernova data and large-scale structure observations all reinforces the evidence for 
dark energy. 

• The simplest model for structure formation, a scale-invariant spectrum of fluctuations, 
is not a good fit to the WMAP data. The WMAP data requires either tensor modes 
or a spectral index with n s < 1 to fit the angular power spectrum. These observations 
match the basic inflationary predictions and are well fit by the predictions of the simple 
m 2 4> 2 model. 

Further WMAP observations and future analyses will test the inflationary paradigm. 
While we do not find convincing evidence for significant non-Gaussianities, an alternative 
model that better fits the low i data would be an exciting development. Within the context 
of the inflationary models, measurements of the spectral index as a function of scale and 
measurements of tensor modes directly will provide a direct probe into the physics of the 
first moments of the big bang. 
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A. SZ Marginalization, Priors, and Likelihood Approximations 

In this Appendix we discuss the affect that various analysis choices have on the estimated 
cosmological parameters, focusing on the six parameter ACDM model. 

The three-year analysis now accounts for the SZ effect by marginalizing over the ampli- 
tude of the SZ contribution, parameterized by the model of Komatsu & Seljak (2002). Specif- 
ically, we evaluate their predicted spectrum using Q m = 0.26, f2 b = 0.044, h = 0.72, n s = 0.97 
and as = 0.8 We define the amplitude of the SZ signal (relative to this model) by A$z and 
marginalize over this parameter with a flat prior, < A$z < 2. This range is based on 
the assumption that the Komatsu & Seljak (2002) (KS) model estimates the true SZ signal 
with an uncertainty of order one. Numerical simulations (Nagai 2005) and analytical studies 
(Reid & Spergel 2006) find a tight correlation between mass and SZ signal, with uncertainties 
that are dominated by the cluster gas fraction. These results support the Komatsu & Seljak 
(2002) approach and suggest that the range of this prior is generous. Afshordi et al. (2005) 
analyze of the SZ signal from 116 nearby clusters in the first-year WMAP data and find that 
the signal from nearby clusters is 30-40% weaker than expected. Since these nearby clusters 
are the dominant source of fluctuations in the WMAP angular power spectrum, this implies 
that A sz < 1. 

We now use the amplitude of the angular power spectrum peak, C220, rather than A as 
the amplitude parameter in the Monte Carlo Markov Chains. This choice of prior leads to 
a slightly lower best fit amplitude. 

Figure 26 shows how the SZ marginalization and the change in amplitude prior affect 
our estimates of cosmological parameters. Except for a s , the effects are all relatively small. 
For as we estimate that roughly half of this change is due to the SZ marginalization, and 
half is due to the change in amplitude prior. Lewis (2006) provides a detailed study of the 
effects of gravitational lensing on the analysis and show that it is smaller than the SZ effect. 

Since the first draft of this paper was circulated, it has been pointed out (Eriksen 
et al. 2006) that the pseudo-Cj-based power spectrum is biased slightly high in the range 
12 < I < 30, which in turn biases the spectral index, n s , slightly low. As a result, we have 
increased the sky map resolution used in our exact TT likelihood module (to N side = 16), 
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which allows us to use this form up to 1=30. We have established that this approach agrees 
well with an independent exact form based on Gibbs sampling methodology (Eriksen et al. 
2006). 

Huffenberger et al. (2006) argue that the point source amplitude quoted in the original 
version of Hinshaw et al. (2006) was too high: they find A ps = 0.011 ± 0.001, while Hinshaw 
et al. (2006) originally quoted A ps = 0.017 ± 0.002. Subsequently Hinshaw et al. (2006) 
have re-analyzed the source contribution and now favor an intermediate value with a more 
conservative error, A ps = 0.014 ±0.003. The original source correction also had the effect of 
biasing the spectral index slightly low. 

Table 15 shows the sensitivity of the spectral index, n s , to the various treatments dis- 
cussed above. The first row (Fiducial) gives the result originally quoted in this paper. The 
second row (No SZ Correction) gives the result if we assume A$z = 0. The third row (No 
Beam Error) shows the effect of assuming the beam window function has no error. The fourth 
row (Lower Point Source Amplitude) shows the change induced by changing the source am- 
plitude (and its uncertainty) from 0.017 ±0.002 to 0.014 ±0.003. The sixth row (Foreground 
Marginalization) shows the effect of disabling marginalization over the foreground emission 
template in the exact likelihood module (up to 1=12). The seventh row (Res 4 Likelihood) 
shows the effect of using the exact likelihood module to 1=30, and the last row shows our 
current best estimate of n s . 

The detailed form of the likelihood function and the treatment of point sources and the 
SZ effect has a ~ 0.5a effect on the best fit slope. We have also used simulations where 
we know the input model to check the effects of including various secondary effects. Lewis 
(2006) provides a detailed study of the effects of gravitational lensing on the analysis and 
show that it is smaller than the SZ effect. 



Table 15: Sensitivity of Slope to Likelihood Model and Prior 



Fiducial (N side = 8 Likelihood, A PS = 0.017) 


0.951 


± 


0.016 


No SZ Correction 


0.954 


± 


0.016 


No Beam Error 


0.953 


± 


0.016 


Reduced Point Source Amplitude (Aps = 0.014) 


0.955 


± 


0.016 


Huffenberger et al. (2006) Point Source Amplitude (A PS = 0.011) 


0.957 


± 


0.016 


Foreground Marginalization 


0.953 


± 


0.016 


N S ide — 16 Likelihood 


0.957 


± 


0.016 


N si de = 16 Likelihood ± Aps = 0.014 ± no SZ correction 


0.960 


± 


0.016 


N side = 16 Likelihood ± A PS = 0.014 


0.958 


± 


0.016 



0.02 0.021 0.022 0.023 0.024 0.08 0.09 0.1 0.11 0.12 0.13 




0.05 0.1 0.15 0.9 0.92 0.94 0.96 0.98 1 




65 70 75 80 85 

H (km/s/Mpc) 



Fig. 26. — The effect of SZ marginalization on the likelihood function. The red 
curve is the likelihood surface for the three-year WMAP data for the power- 
law ACDM model with A S z = 0. The black curve is the likelihood surface after 
marginalizing over the amplitude of the SZ contribution. 
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B. Trispectrum methodology 

B.l. Predicted Trispectrum Signal 

We consider here a model that generates a non-trivial trispectrum, but no bispectrum 
signal. We assume that the gravitational potential $, is a product of two independent 
Gaussian fields, and tp: 

$(£)= <f>(x) [l + g NL ^(x)\ (Bl) 
where g^L characterizes the strength of the non-linear term. 

Following Komatsu et al. (2005) approach for the bispectrum, extended recently to the 
trispectrum in Kogo & Komatsu (2006), the observed temperature multipoles are: 

».-k/-H*Wr)+^ (B2) 
where bg is the beam, n^ m is the noise and o^(r) is the radiation transfer function: 

a e (r) = ^J k 2 dkg Tl {k) 3 t{kr) (B3) 

The non-linear coupling term generates a second order term: 

aim = n im + b e J r 2 dra t {r) <p im (r) + <pi> m ir)4} illmll (r)C e e '™' e '' m ''j (B4) 

where 



(2£+ l)(2f + l)(2f + 1) V / \ m m' m" 

This term does not have any effect on the bispectrum as < 3 >= and < ^ >= 0. 
However, it does have a non-trivial effect on the trispectrum. 

As with gravitational lensing (see Hu (2001)), the largest trispectrum term is the diag- 
onal term, T//(0) =< CiPt > — 3 < Ce > 2 ■ This term would generate an excess in the x 2 of 
the fit of the model to the data: 



T'M = 9 2 NL C e J r 2 dr J f 2 dfa e (r)a e (f) 



ml'm! l"m" H' m ' l"m" 
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We can then use 

< ^'m'( r )^ w (0 >= S e r e ,S m >m> J k 2 dkP,p(k)je(kr)j e (kr) (B7) 
and the equivalent relationship for ip to rewrite the trispectrum as 

Tg(0) = 9 2 NL ^Cj (B8) 

where 



(II 1)11 " 2 



r 2 dr / r 2 drae(r)ae(r)j e >(kr)ji»(k f r)ji>(kr)ji»(k'r). (B9) 





While this full integral is numerically intractable, we approximate the surface of last scatter 
as a thin screen so that 

a im = bt$em(r*)a e + n lm (BIO) 
then, the trispectrum coupling term reduces to 

2 

6 = ^r- / k 2 dkP^(k)j 2 (kn) J k' 2 dk'P4k')fAk'n) ( J £ ^ ) . (Bll) 
Recall that in this limit, 

C< = a, 2 y k 2 dkP^k)jj{kr,) (B12) 

Thus, 

6 = ^J^W^'r.) ( J J £)' (B13) 

The amplitude of ^£ is, thus, roughly the variance in the ip field on the scale r*jl. Note that 
£i is a positive definite quantity so that T//(0) should be always positive. 



B.2. Detecting the Non-Gaussian Signal 

If we assume that ^ is constant, then we can follow Hu (2001) and compute an optimal 
quadratic statistic. We approximate the optimal statistic as ^Z^T/^T? — N 2 ) 2 , where 
is a smoothed map {e.g. a res 7 map) and we use the approximation that Cg = A/ 1(1 + 1). 
This has the advantage that we can easily compute it and has well-defined noise properties. 
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B.2.1. Practical implementations 
We define for this purpose the dimensionless Q statistic as 

Q = ^ W p,b i r ^P,bi W P,b2^P,b2 W P,b-A,bs W P,bAM 

P,&1,&2,&3,&4 

bi,b 2 ,b 3 ,b4 Pi P2 

(B14) 

where fej refers to various bands (Q, V, and W for yrl, yr2 and yr3) that are all distinct for 
a single term so that the noise bias is null for this statistic, w Pjbi is a particular pixel weight 
(we will consider it equal unity first) and T b is a filtered map defined as 

f pb = (B15) 

\JY,q T qb 

T pb = Y^f^fmbY^p) (B16) 

tm 

where a^'s are the spherical harmonic coefficients of the masked sky. We use at this level the 
Kpl2 mask to hide the brighter part of the galaxy (and potentially the brighter point sources) 
and ignore cut sky effects in considering those pseudo-a^s. But when computing the sum 
over pixels in Q, we consider only pixels outside the Kp2 area. The obvious advantage of 
this simple real space statistic is its ability to handle inhomogeneous noise and to localize its 
various contributions in real space. The second term in the definition of Q aims at subtracting 
off the Gaussian unconnected part, so that if the T fields are homogeneous Gaussian fields, 
we obtain {Q ps ) = 0. 

The exact nature of f e will depend on the source of the signal. For example, point 
sources do contribute to all n-points functions in real or harmonic space and are as such 
visible in the power spectrum, bispectrum and trispectrum. The first two have been used to 
set limits and corrections. 

Should we want to isolate the point sources contribution with the Q statistic, we would 
proceed in the following manner. The point sources power spectrum is well approximated 
by a constant, white noise like, power spectrum C v t s (see Komatsu et al. (2003)). Given the 
measured power spectrum, C e = C e bj + C^ s bj + N e , the Wiener like filter to reconstruct point 
sources is to ff s = b\jC t Note that this filter would be optimal only if point sources were 
drawn from a Gaussian distribution, which is not true. We can however expect it to be close 
to optimal. 
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In order to constrain the CMB contribution to the trispectrum and constrain g NL close 
to optimality, we will set f n and f e2 to the Wiener filter for the CMB field, f £1 = f e3 = 
c theory b 2j C measured and = = ^ + i)/ fl> where Cf eory is the best fit model angular 
power spectrum and C^ ieasured is the measured raw power spectrum including the signal and 
the noise and not corrected for the beam window function. 

We restrict ourselves to a unit weighting which is nearly optimal in the signal dominated 
regime where we draw our conclusions from, i.e. at resolution lower than N side = 64 



B.2.2. Explicit relation to the trispectrum 
Ignoring the weights, it is easy to show using the relation recalled in the next section 

that 

(T 1 (n)T 2 (n)T 3 (n)T 4 (n)) c = J2 T ^p T ^p ( B17 ) 

v 

Kn^h^imjlmjcf dtt(fl)Y iin ^ 



l\ mi <2m2 £3013^4 "M 



(B18) 



l (2f 1 + l)(24 + l) / (2l 3 + l)(2l4 + l) r£0 r L0 r LM r LM 



Air(2L + 1) V 47r(2L + 1) h0(-20 ^O^O (-im 1 e 2 m 2 ^tzmzlAm^, 

X (^imi^2m2^3m3^4m4/c ■ (B19) 



It is then easy to relate to standard expression for the connected part of the trispectrum as 
in Hu (2001) and Komatsu (2001). 
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